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ABSTRACT 

The redshifts of all cosmologically distant sources are expected to experience a small, sys- 
tematic drift as a function of time due to the evolution of the Universe's expansion rate. A 
measurement of this effect would represent a direct and entirely model-independent deter- 
mination of the expansion history of the Universe over a redshift range that is inaccessible 
to other methods. Here we investigate the impact of the next generation of Extremely Large 
Telescopes on the feasibility of detecting and characterising the cosmological redshift drift. 
We consider the Lyman a forest in the redshift range 2 < z < 5 and other absorption lines in 
the spectra of high redshift QSOs as the most suitable targets for a redshift drift experiment. 
Assuming photon-noise limited observations and using extensive Monte Carlo simulations 
we determine the accuracy to which the redshift drift can be measured from the Lya forest 
as a function of signal-to-noise and redshift. Based on this relation and using the brightness 
and redshift distributions of known QSOs we find that a 42-m telescope is capable of unam- 
biguously detecting the redshift drift over a period of ^20 yr using 4000 h of observing time. 
Such an experiment would provide independent evidence for the existence of dark energy 
without assuming spatial flatness, using any other cosmological constraints or making any 
other astrophysical assumption. 

Key words: cosmology: miscellaneous - intergalactic medium - quasars: absorption lines. 



1 INTRODUCTION 

The universal expansion was the first observational evidence that 
gener al relativ it y mig ht be applicable to the Universe as a whole. 
Since lHubbld 's ( 1 1929" ) discovery much effort has been invested into 
completing the basic picture of relativistic cosmology. The central 
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question is: what is the stress-energy tensor of the Universe? As- 
suming homogeneity and isotropy reduces this question to: what is 
the mean density and equation of state of each mass-energy com- 
ponent of the Universe? Since these parameters determine both 
the evolution with time and the geometry of the metric that solves 
the Einstein equation, one can use a measurement of either to in- 
fer their values. Over the past decade the successes on this front 
have reached their (temporary) culmination: observations of the 
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Cosmic Microwave Backgr ound (CMB 



type la supernovae (SNIa; iRiess et al 



Spergel et alj2003 



2004 



Astier et al 



the large-scale galaxy distribution ('Peaco ck et alJl200ll : ICo 



2007), 



200€), 



e et 



[2005; Eisenstein et al. 2005) and others now provide answers of 
such convincing consistency and accu racy that the ter m 'precision 
cosmology' is now commonplace (e.g. |Primackl200l) . 

By far the most unexpected result of this campaign was the 
discovery that the expansion of the Universe has re cently begun 
accelerating tRiess et ai] |l998l :l Perlmutter et iD[l999|). The phys- 
ical reason for this acceleration is entirely unclear at present. 
Within relativistic cosmology it can be accommodated by mod- 
ifying the stress-energy tensor to include a new component with 
negative pressure. In its simplest incarnati on this so-called dark en- 
ergy is the cosmological constant A (e.g. ICarroU. Press & Turnetl 
Il992h . i.e. a smooth, non-varying component with equation of 
state paramete r w = — 1. It could also be time variable (e.g. 
[Overduin & CooDerstockl ll998), either b ec ause o f a specific equa- 
tion of state (x CDM, e.g. rTumer & White! 1 19971 : phantom energy, 
ICaldwellll2002l) or because the componen t consists of a dynami - 
cal scalar field evolving in a potential (e.g. lRatra & Peebleslll988h . 
More generally it may be an inhomogeneous, time varying com- 
ponent with —^'^w<Q (where w may also vary with time), 
some times called quintessence (e.g. ICaldwell. Dave & Steinhard3 
Il998l) . or even a c omponent with an exotic equation of stat e 
(e.g. Chaplygin gas. iKamenshchik. Moschella & Pasauiej|200l[) . 
Alternatively, instead of modifying the stress-energy tensor 
one can also modify gravity itself to explain the acceleration 
(e.g. iDeffavet, Dvali & Gabadadzd I2OO2I : Ipreese & Lewij I2OO2I : 
ICarroU et alj|2005l) . Although all curr ent observations are consis- 
tent w ith a cosmological constant (e.g. lSeliak. Slosar & McDonald! 
2006h many more un familiar models are not ruled out (e.g. 



Capozziello et al.!!2005!) . 



Probably the best way to probe the nature of the accelera- 
tion is to determine the expansion history of the Universe jUindeJ 
[2003; Seo & Eisenstein 2003). Observables that depend on the ex- 
pansion histor y include distances and the linear gro wth of density 
perturbations dLinder & Jenkinsl 2003 : LindeJ ioosi). and so SNIa 
surveys, weak lensing 1 HeavensT 2003 : Jain & Taylor! !2003l) and 
bary on acoustic oscillations in the galaxy power spectrum (BAO; 
!Seo & Eisenstein 2003; Wang 2006) are generally considered to be 
excellent probes of the acceleration. 

In practice, however, extracting information on the expansion 
history from weak lensing and BAO requires a prior on the spa- 
tial curvature, a detailed understanding of the linear growth of den- 
sity perturbations and hence a specific cosmological model. Given 
the uncertain state of affairs regarding the source of the accelera- 
tion, and given that even simple parameterisatio ns of dark energ; 
properties can resu lt in misleading conclusions jMaor et 



K energy 
ai] !200li ; 



!Bassett et alj2004) . these are conceptually undesirable features and 
several authors have pointed out the importance of taking a cosmo- 
graphic, mo del-independent approach to determining the expansion 
histor y (e.g. !Wang & Tegmarkl!2005! ; !johnl!2005l ; !ShaDiro & Turner! 
!2006l) . Using SNIa to measure luminosity distances as a function 
of redshift is conceptually the simplest experiment and hence ap- 
pears to be the most useful in this respect. The caveats are that dis- 
tance is 'only' related to the expansion history through an integral 
over redshift and that one still requires a prior on spatial curvature 
JCaldwell & Kamionkowski 2004). 

Here we will revisit a method to directly measur e the expan- 
sion history that was first explored bv !Sandag3 l ! 19621 with an ap- 
pendix by McVittie). He showed that the evolution of the Hubble 
expansion causes the redshifts of distant objects partaking in the 



Hubble flow to change slowly with time. Just as the redshift, z, is 
in itself evidence of the expansion, so is the change in redshift, i, 
evidence of its de- or acceleration between the epoch z and today. 
This implies that the expansion history can be determined, at least 
in principle, by means of a straightforward spectroscopic monitor- 
ing campaign. 

The redshift drift is a direct, entirely model-independent mea- 
surement of the expansion history of the Universe which does not 
require any cosmological assumptions or priors whatsoever. How- 
ever, the most unique feature of this experiment is that it directly 
probes the global dynamics of the metric. All other existing cos- 
mological observations, including those of the CMB, SNIa, weak 
lensing and BAO, are essentially geometric in nature in the sense 
that they map out space, its curvature and its evolution. Many of 
these experiments also probe the dynamics of localised density per- 
turbations but none actually measure the global dynamics. In this 
sense the redshift drift method is qualitatively different from all 
other cosmological observations, offering a truly independent and 
unique approach to the exploration of the expansion history of the 
Universe. 

Following the original study by ISandage! Il962l) . the redshift 
drift and its relevance to observ ational cosmo lo gy were also dis 



cussed byJMcVittie (1965), W einberg! (!l972 ). lEbert & Triimper 
1975h . iRUdigen (|198Q), iPeacock! il999!). kakamura & Chiba 



I999I) , !Zhu & Fuiimotol 1 2004). !Corasaniti. Hute rer & Melchiorri 



1 2007) an d' Lake! ( !2007hlLakel ( |198I|) gave equations expressing the 
deceleration and matter density parameters, go and Q.m, in terms of 
z, z and z. An excellent expose of t he equations relevant to red- 
shift evolution was also presented by !Gudmundsson & Biornsson! 
(2002). In addition, these authors investigated the redshift drift in 
the presence of quintessence, w hile other non-standa r d dark en- 
ergy models were considered by !Balbi & Ouercellinil ( !2007h and 
!zhang et al.! ( !2007|k The case of Dicke-Brans-Jordan cosm ologies 
was scrutinised by lRudigei!fl982b and Partovi & Mashhoon' fl984 
studied the effects of inhomogeneities. Seto & Cooray (200^ sug- 
gested that a measurement of the anisotropy of z could be used to 
constrain the very low-frequency gravitational wave background. 
Several of these authors have pointed out the superior redshift ac- 
curacy achieved in the radio regime compared to the optical, and 
[Oavis & May ( 1978) entertained the possibility of using a 21 cm 
absorption line at z — 0.69 in the radio spectrum of 3C 286 for a 
z measurement. The detrimental effects of peculiar velocities and 
ac celerations, which may s wamp the cosmic signal, were disc ussed 
by Phiiiippi ( !l982h . !Lakel ( !l982l) and !Teubei! jl986!) . Finally. !Loebl 
^1998) first proposed the Lyman a (Lya) forest as an appropri- 
ate target for a i measurement (acknowledging D. Sasselov for the 
idea) and assessed the prospects for a successful detection in the 
context of currently existing observational technology. All except 
the last of these studies concluded that a i measurement was be- 
yond the observational capabilities of the time. 

The purpose of the present paper is to examine the impact 
of the next generation of 30-60-m Extremely Large Telescopes 
(ELTs) on the feasibility of determining z. The key issue is ob- 
viously the accuracy to which one can determine redshifts. In the 
absence of systematic instrumental effects this accuracy depends 
only on the intrinsic sharpness of the spectral features used, the 
number of features available and the signal-to-noise ratio (S/N) at 
which they are recorded (assuming that the features are resolved). 
In a photon-noise limited experiment the latter in turn depends only 
on the flux density of the source(s), the size of the telescope, the 
total combined telescope/instrument efficiency and the integration 
time. In this paper we will investigate this parameter space in detail. 
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expanding on previous work in several ways: in Section[3]we con- 
firm the usefulness of the Lya forest by quantifying the peculiar 
motions of the absorbing gas using hydrodynamic simulations of 
the intergalactic medium (IGM). In Section|4]we use Monte Carlo 
simulations of the Lya forest to quantify how its properties trans- 
late to a radial velocity accuracy, and we consider the usefulness of 
other absorption lines in Section |5] These results are then used in 
Section|7]where we explore the observational parameter space and 
realistically assess the feasibility of a i experiment with an ELT. 
Finally, we summarise our findings in Section[8] 



2 MEASURING THE DYNAMICS 

In any metric theory of gravity one is led to a very specific form of 
the metric by simply assuming that the Universe is homogeneous 
and isotropic. The evolution in time of this so-called Robertson- 
Walker metric is entirely specified by its global scale factor, a{t). 
The goal is to measure or reconstruct this function. Recall that we 
observe the change of a with time by its wavelength-stretching ef- 
fect on photons traversing the Universe. A photon emitted by some 
object at comoving distance x at t™e ton and observed by us at 
tobs suffers a redshift of 



Big Bang 



1 + z{tohB,t cm / 



a{tob 



(1) 



(where, obviously, only two of the three variables x, iem and tohs 
are independent). If it were possible to measure not only a pho- 
ton's redshift but also its tcm, then the problem would be solved by 
simply mapping out the present-day relation between redshift and 
look-back time, i.e. z{tobs ~ to, tcm), where to denotes today. 

A different approach to the problem is to consider how the 
redshift of an object at a fixed comoving distance x evolves with 
fobs, i.e. to consider the function 2|^(fobs). For a given x- tohs 
determines fom and so we have dropped the dependence on ton. In 
principle, it is possible to map out (fobs) (at least for fobs > to), 
and this would be the most direct determination of a{t), no matter 
which object is used. However, a full characterisation of Z|;^(fobs) 
would require observations over several Gyr. Over a much shorter 
timescale, Afobs, one can at most hope to determine the first order 
term of the Taylor expansion 



dz|- 
dfoi 



-{to 



•2|^(^obs H" ^^obs 



2|X 



(^obs 



Afo 



(2) 



As we will see presently, it turns out that measuring d^i^/dfobs 
is in fact sufficient to reach our goal of reconstructing a{t). By 
differentiating equation (TJ with respect to fobs we find 



■(fobs) = [1 + Z|x(tobs)]-H'(fobs) - -H'(fcm), 



(3) 



dfobi 

where we have used that dfobs = [1 + z(fobs)] dfcm for a fixed x, 
and H = aa~^ . Evaluating at fobs = fo, replacing the unknown 
fern with its corresponding redshift, and dropping the reminder that 
we are conside ring the redshif t of an object at a fixed distance x we 
simply obtain jMcVittiel 19621) : 



dz 

dfobs 



(fo) = {\ + z)Ho~a{z). 



(4) 



i is a small, systematic drift as a function of time in the redshift 
of a cosmologically distant source as observed by us today. This 
effect is induced by the de- or acceleration of the expansion, i.e. 
by the change of the Hubble parameter H. Since Ho is known 
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Figure 1. Redshift evolution of three objects with present-day redshifts 
z(to) = 0.5, 3 and 8 as a function of time of observation, for three dif- 
ferent combinations of Qm and as indicated. For each case the dotted 
lines indicate the Big Bang, fobs is shown relative to the present day, fo. 



jFreedman et aljfaoil) . this drift is a direct measure of the expan- 
sion velocity at redshift z. Measuring i for a number of objects 
at different z hence gives us the function a{z). The point is that 
given a{z) and d{z), one can reconstruct a(f). A measurement of 
z{z) therefore amounts to a purely dynamical reconstruction of the 
expansion history of the Universe. 

Predicting the redshift drift z{z) requires a theory of gravity. 
Inserting the Robertson-Walker metric into the theory's field equa- 
tion results in the Friedman equation, which specifically links the 
expansion history with the densities, fii, and equation of state pa- 
rameters, Wi, of the various mass-energy components of the Uni- 
verse. In the case of general relativity the Friedman equation is 
given by: 



H{z) = Ho 



J2n,{i+z: 



+ nk{i + zy 



(5) 



where fifc = 1 — ^ Here we consider only two components: 
cold (dark) matter (CDM) with wm ~ and dark energy in the 
form of a cosmological constant with wa ~ —1. In Fig.[T]we plot 
2|x(iobs) for three different objects [chosen to have z{to) = 0.5, 3 
and 8] and for three different combinations of Q,m and f^A, where 
we have also assumed Ho — 70 ftro km s^^ Mpc^^. For each 
object, the redshift goes to infinity at some time in the past when 
the object first entered our particle horizon. If SIa = the redshift 
continually decreases thereafter as the expansion is progressively 
slowed down by flu - Hence, in this case the redshifts of all objects 
are decreasing at the present time. 

However, if Qa the initial decrease is followed by a 
subsequent rise due to A relieving matter as the dominant mass- 
energy component and causing the expansion to accelerate. The 
turn-around point may lie either in the past or in the future, de- 
pending on the object's distance from us; i.e. if SIa 7^ an object's 
redshift may be either increasing or decreasing at the present time. 
For distant (nearby) objects, the Universe was mostly matter (A)- 
dominated during the interval [fcm, fo] and hence underwent a net 
decelerat ion (acceleration), resulting in i < (>) at the present 
time (see lGudmundsson & Bi6mssoi]||2002l for a more detailed dis- 
cussion). 
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Figure 2. The solid (dotted) lines and left (right) axis show the redshift drift 
z (i)) as a function of redshift for various combinations of Hm and Ha as 
indicated. The dashed line shows i for the case of dark energy having a 



constant wx 



■ (and Hm, fix = 0.3,0.7). 



These features are evident in Fig. (2] where we plot the ex- 
pected present-day redshift drift, z{z), for various values of SIm 
and Q,A (solid lines), and for a case where the dark energy «; 7^ — 1 
(dashed line). The redshift drift is also shown in velocity units (dot- 
ted lines), where v — c z {1 + z)~^ . As, noted above, the existence 
of a redshift region where i > is the hallmark of Q.^ 7^ 0. Note 
also the scale of Fig.|2] At 2 = 4 the redshift drift is of order 10"' 
or 6 cm s"^ per decade. For comparison, the long-term accuracy 
achieved in extra-solar planet searches with the high-resolution 
echelle spectrogra ph HARPS on the ESO 3.6-m Telescope is of 
order 1 m s"^ (e.g. lLovis etal]|2005h . 



3 CHOOSING AN ACCELEROMETER 

A priori, it is not at all obvious which spectral features of which 
set or class of objects might be best suited for a z measurement. 
Clearly though, potential candidate targets should boast as many 
of the following desirable attributes as possible, (i) They should 
faithfully trace the Hubble flow. Although peculiar motions are 
expected to be random with respect to the Hubble flow, the addi- 
tional nois e introduced by them could p otentially conceal the cos- 
mic signal ("Phillipp sl 198a : lTeubeil 198 9; but see also Appendix lAt. 
(ii) The targets should have the sharpest possible spectral features 
to minimise the error on individual redshift measurements, (iii) The 
number of useful spectral features per target should be as high as 
possible in order to maximise the amount of relevant information 
per unit observing time, (iv) The targets should be as bright as pos- 
sible and finally, (v) they should exist over a wide redshift range, 
and particularly at high z, where the signal is expected to be largest. 

Clearly, several of these features are in conflict with each 
other. Requirements (i) and (ii) are conflicting because sharp spec- 
tral features require cold material which is generally found in dense 
regions inside deep potential wells, which in turn generate large 
peculiar accelerations. Similarly, point (i) clashes with point (iv) 
because high intrinsic luminosities require very massive and/or 
highly energetic processes, again implying deep potential wells. 
Since very dense regions are relatively rare there is likewise ten- 



sion between requirements (ii) and (iii). Finally, the demands for 
brightness and high redshift are also difficult to meet simultane- 
ously. 

Hence, it seems that it is impossible to choose a class or set 
of objects that is the ideal target, in the sense that it is superior to 
every other class or set in each of the above categories. However, 
there is one class of 'objects' that meets all but one of the criteria. 



3.1 The Lyman a forest 

The term 'Lya forest' refers to the plethora of absorption lines 
observed in the spectra of all quasi-stellar objects (QSOs) short- 
wards of the Lya emission lineQ Almost all of this absorption 
arises in intervening intergalactic H I between us and the QSO (see 
[Ranch 1998 for a review). Since the absorbing gas is physically 
unconnected with the background source against which it is ob- 
served we elegantly avoid the conflict between requirements (i) 
and (iv) above. However, as the gas is in photoionization equi- 
librium with the intergalactic ultrav iolet (UV) background, its 
temperature is of o rder 10* K (Theun s. Schave & Haehnelt 200(3; 
ISchave et alj2000h . Consequently, the absorption lines are not par- 
ticularly sharp and the typical li ne width is ~30 km s~^ (e.g. 
iKim. Cristiani & D'Odoricoll200ll) . On the other hand, QSOs are 
among the brightest sou rces in the Universe and exist at all red- 
shifts out to at least ~6 jFan et al . 2006). Furthermore, each QSO 
spectrum at z > 2 shows on the order of 10^ absorption features. 

In the following we will consider the question whether, apart 
from the Hubble expansion, other evolutionary processes acting on 
the absorbing gas might also significantly affect the measured po- 
sitions of the absorption lines over the timescale of a decade or 
so. First and foremost is the issue to what extent the Lya forest is 
subject to peculiar motions. 



3.1.1 Peculiar motions 



Cosmological hydrodynam ic sim ulat ions 1 Cen et al] 
Zhang. Arminos & Normar 1 19951; iHernguist et al! 



elling (e.g 
observations 



1994 



1996 



1999 ; ICen & Simcod Il997l ; Icharlton et al.l 



Miralda-Escude et al.l 

19971 ; iTheuns et al.l Il998l; Izhang et alj 1 19981), an alytic mod- 



Bi&DavidserJ 1 19971 : 



IViel et al.l l2002h and 

o f the sizes and shapes of the absorbers 

teechtold et ai] Il994l; [Pinshaw et al.l 1 19941; ISmette et alj 1 19951; 
Charlton. Churchill & Linded 1 1 9951; [P'Odorico et alj 1 19981; 
Liske et aL; 2000i; Rol l ijide et alj 120031; (secker. Sargent & RauchI 
2004; ICoppolani et alj|200i; see also iRauch & Haehnellll 19951) all 
suggest that the Lya forest absorption occurs in the large-scale, 
filamentary or sheet-like structures that form the cosmic web. 
These are at most mildly overdense and are participating, at least 
in an average sense, in the general Hubble expansion. However, 
at least some fraction of the absorbing structures (depending on 
redshift) must be expected to have broken away from the Hubble 
flow and to have begun collapsing under the influence of local 

gravitational potential wells. 

This exact issue was investigated bv lRauchet"ai] ( l2005b who 
studied the distribution of velocity shear between high-redshift Lya 
absorption common to adjacent lines of sight separated by 1 to 



^ Since the Lya forest is only observable from the ground for 2 > 1.7 
we only consider the high-z Lya forest in this paper. We also exclude the 
higher column density Ly limit and damped Lya (DLA) absorbers from the 
discussion. 
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300 h^Q kpc. They showed that the observed shear distributions 
were indeed in good agreement with the absorbing structures un- 
dergoing large-scale motions dominated by the Hubble flow. In 
fact, the distributions could be reproduced very convincingly by ar- 
tificial pairs of spectra created from a hydrodynamical simulation 
of the IGM. This result instills us with further confidence that such 
simulations accurately capture the kinematics of the gas responsi- 
ble for the Lya forest. 

Thus assured, we will now explicitly examine the peculiar 
velocities and accelerations of the absorbing gas in a hydrody- 
namic simulation of the IGM pr oduced by the parallel TreeSPH 
code GADGET-2 ( ISpringelir2005l) . We have used this code in its 
TreePM mode in order to speed up the calculation of long-range 
gravitational forces. The simulations were performed with periodic 
boundary conditions and with 400^^ dark matter and 400^ gas par- 
ticles in a box of 60 /i^qq Mpc comoving size. Radiative cooling 
and heating p rocesses were followed using an imp lementation sim- 
ilar to that of lKatz. Weinberg & Hemguistl ( Il99q) for a primordial 
mix of hydrogen and he lium. The UV background was taken from 
lHaardt&Madaul l ll996h . In order to further increase the speed of 
the simulation we applied a simplified star formation criterion: all 
of the gas at overdensities > 10"^ times the mean density and with 
temperature < 10^ K was turned into stars. The cosmological pa- 
rameters were set to = 0.26, f^A = 0.74, f^b = 0.0463, 
Us = 0.95, (78 = 0.85 and Ho = 72 km s"^ Mpc"\ where 
flb is the baryonic density parameter, and and erg are the spec- 
tral index and amplitude of the linear dark matter power spectrum, 
respectively. These values are in excellent agreement with recent 
joint analyses of CM B , SNIa, galaxy clus t ering and hya forest data 
dSpergel et alj |2007| ; ISeliaketalJ I2OO6I : Iviel. Haehnelt & LewisI 
I2OO6I). TheACD M transfer function was computed with CMBFAST 
jSeliak & Zaldarriaga 1996,). The above cor r espon ds to the simula- 
tion series B2 of Viel. Haehnelt & Springell ( l2004h which has been 
widely used for cosmological studies. Moreover, this is the exact 
same simulation that su ccessfully reproduc ed the velocity shear 
distributions observed bv lRauch et al.l l2005h . 

We pierced our simulation box with 1000 random lines of 
sight (LOS) and noted the physical properties of the absorbing 
gas along these lines at three different redshifts (2 = 2, 3 and 
4). In Figs. [3] and |4] we plot the distributions of peculiar velocity 
and acceleration, respectively. The acceleration of the gas is com- 
puted from its velocity by dividing th e latter by a dynamical time 
tdyn = (Gp)"^''^ (e.g. lSchav3l200ll) . where p is the total matter 
density. Evidently, the distributions do not evolve rapidly with red- 
shift. Since a given parcel of gas cannot be decelerated, peculiar 
velocities cannot decrease with time and hence their distribution 
shifts to slightly larger values at lower redshifts. The average pecu- 
liar acceleration, on the other hand, decreases with time because by 
volume the Universe is dominated by low-density regions in which 
the gas density decreases as the Universe expands. 

Peculiar motions are expected to be randomly oriented with 
respect to the LOS. Hence, when averaging over a large num- 
ber of individual i measurements, peculiar motions will only in- 
troduce an additional random noise component but no systematic 
bias. In Appendix lAl we explicitly derive an expression for the ob- 
served redshift drift in the presence of peculiar motions. This ex- 
pression can be used to translate the distributions in Figs. [3] and |4] 
into the corresponding error distribution on 2. We find that for a 
decade-long experiment the error due to peculiar motion is of or- 
der ~10~^ cm s~^. This must be compared to the error induced by 
photon noise for an individual z measurement from a single absorp- 
tion line. Clearly, if we wish to be able to detect the redshift drift 
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Figure 3. Probability distribution functions (PDFs) of the peculiar velocity 
of the absorbing gas along 1000 random lines of sight through our simula- 
tion box at three different redshifts as indicated. Note that we are not using 
the value of the velocity's component along the LOS, but rather the modulus 
of the full 3-dimensional velocity. 
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Figure 4. As Fig.[3]for the peculiar acceleration. 



over a decade or so the overall accuracy of the whole experiment 
has to be of order ^ 1 cm s~ ^ . Since this will be achieved using hun- 
dreds of absorption lines the error on an individual 2 measurement 
will be at least a factor of ~10 larger. The error due to peculiar mo- 
tions is therefore sub-dominant by at least ~4 orders of magnitude. 
Hence we conclude that peculiar motions will have no detrimental 
impact whatsoever on a redshift drift experiment targeting the Lya 
forest. 
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3.1.2 Galactic feedback 

So far we have only considered peculiar motions induced by grav- 
ity. In principle, there are of course a number of non-gravitational 
ways of imparting kinetic energy to the absorbing gas, mostly 
involving galactic feedback. The main reason to believe that 
galactic feedback must have had a far-reaching impact on the 
Lya forest is its early and widespread low-l evel metal enrich- 
ment, even at fairly low column densities (e 



Ellison et al." 2000: 'Schave et alj |2000l 120031: 1 Aguirre et al 



Cowie et al 



1995 



2004 



Simcoe, Sargent & Rauch 2004)). Although very low-level enrich 
ment could be achieved in situ with Population III stars (e.g. 
iGnedin & Ostrikeilll997l) the general consensus is that some sort 
of mechanism is required to transport metals from (proto-)galaxies 
into the IGM. In contrast, no consensus has been reached as to 
which of the possible mechanisms (or combination thereof) is the 
correct one. Candidates are mergers and tidal interactions (e.g. 
lGnedinl ll998'), ram pressure stripping, radiation pressure on dust 
grains (e^. Aguirre et al. 2001), photo-evaporation during reion- 
ization ( IBarkana & Loeblll999l) and SN-driven winds, either from 
low-mass, (pre-)galactic halos at z ~ 10 (e.g. iMadau et al.1l200lh 
or from massive starbursting galaxies at z < 5 (e.g. I Aguirre et aU 
l200lUAdelberger et alj2003l) . From our point of view the last pos- 
sibility is the most worrying as it has the highest potential of sig- 
nificantly altering the kinematic structure of the IGM at the time of 
observation. 

There exists persuasive evidence of the existenc e of galactic 
superwinds from Lyman break galaxies at 2: ~ 3 (e.g. lPettini et al.l 
[200 1) and it is likely that some fraction of stro ng metal absorp- 
tion l i nes are connected w ith these structures (e.g. lAdelberger et al.l 
I2OO3I : [ Simcoe et al.ll200^ . On the other hand, there is no obser- 
vational evidence at all that superwinds are significantly stirring 
up the high redshift Lya forest at the time we observe it. Using 
column density a nd optical d e pth di fferences across a close pair 
of lines of sight iRauch et alj hoOlh found no indication in the 
forest's small-scale density structure for widespread recent distur- 
bances. Similarly, as discussed above, the velocity shear between 
adjacent lines of sight is entirely explained by the Hubble flow and 
gravitational instability jRauch et al. 2005), leaving little room for 
non-gravitationally induced motion. Indeed, the simple fact that 
the aforementioned hydrodynamical simulations - which did not 
include any galactic feedback - were so successful in reproduc- 
ing the observed properties of the Lya forest, including its line 
broadening distribution and clustering, raises the question of how 
significant amounts of feedback could be integrated without up- 
setting the existing agreeme nt between the models and the data 
llh euns. Mo & Schav j|200lh . It seems that the volu me filling fac 



tor of galactic superwinds is limited to a few per cent jTheuns et al 
2OO2I; iDesiacQues et alj |2004| : iPieri & Haehneij |2004| ; ICen etal 



20051; iBertone et al.ll2005l) . We conclude that, whatever the process 



of metal enrichment may be, there is currently no reason to believe 
that it has a wholesale effect on the kinematics of the general IGM 
(as probed by the Lya forest) at the time of observation. 



3.1.3 Optical depth variations 

Consider the gas responsible for a given Lya forest absorption fea- 
ture. If the physical properties of the gas change over the timescale 
of a decade or so this will cause a variation of the feature's optical 
depth profile. Potentially, this could in turn lead to a small shift in 
the feature's measured position and hence mimic a redshift drift. 
The precise magnitude of this additional i error component will 



depend on the method used to extract the signal but we can gain 
an impression of the relevance of the effect by comparing the ex- 
pected optical depth variation to the apparent optical depth change 
(at a fixed spectral position) induced by the redshift drift. 

The gas properties we consider here are density, temperature 
and ionizati on fraction. In the fluctuating Gunn Peterson approxi- 
mation (e.g. iHui. Gnedin & Zhang|ll997 ) the Lya optical depth is 
related to these quantities by r oc (1 + 5)^ T^"'^ F^^, where S 
and T are the gas overdensity and temperature, respectively, and F 
is the photoionization rate. The interplay between photoionization 
heating of the gas and adiabatic cooling leads to a tight relation be- 
tween temperatur e and density, which c an be well approximated by 
r = ro(l + Sy ( lHui& Gnedi"i]|l997h . Hence we obtain 



r oc (1 + (5) 



2-0.77 



Tn 



(6) 



How do these quantities evolve with time? According to linear the- 
ory density perturbations grow as ( 1 + z ) ~ ^ . Although exactly true 
only for an Einstein-de Sitter Universe, this represents an upper 
limit in open and flat models with a cosmological constant. Hence 
we will err on the side of caution by adopting this growth factor in 
the following. The evolution of the temperature-density relation can 
be gleaned from figure 6 of Schave et al. (2000): dj/dz = —0.1 
and dTo/dz = 0.25 To (where we have conservatively used the 
steep evolution betw een redshifts 2 and 3). We take the evolution 
of F from figure 7 of iBolton et al] ( l2005h : dF/dz = -0.3F. Tak- 
ing the appropriate derivatives of equation inserting the above 
values and translating a decade in our reference frame to a red- 
shift difference dz at 2: = 3 we find an optical depth variation of 
dr = 3 X 10"^ /r, where / « 0.2, 0.2, 0.05, 0.3 for the variation 
of 5, To, 7 and F, respectively. As we will see in Section l4!4l these 
values for dr are at least ~2 orders of magnitude smaller than the 
typical optical depth changes due to the redshift drift. We therefore 
conclude that changes in the physical properties of the absorbing 
gas are not expected to interfere with a z measurement from the 
Lya forest. 

3.2 Molecular absorption lines 

Before we move on to investigate the details of a redshift drift ex- 
periment using ELT observations of the Lya forest, let us briefly di- 
gress here to consider a very different z experiment using another 
future facility. We have explored in some detail the possibility of 
using the Atacama Large Millimeter Array (ALMA) to measure 
z from rotational molecular transitions seen in absorption against 
background continuum sources. Going to the (sub-)mm regime has 
the advantages of potentially very high resolution and less photon 
noise for a given energy flux. Furthermore, molecular absorption 
lines can be very sharp, with line widths as low as < 1 km s^^. 
However, the molecular gas in nearby galaxies is strongly concen- 
trated towards the central regions. Hence we must expect the gas to 
be subject to peculiar accelerations similar in magnitude to the cos- 
mological signal. That in itself would not necessarily be problem- 
atic as long as we had many individual 2 measurements from differ- 
ent objects over which to average. Unfortunately, it seems unlikely 
that this will be the case. At present, rotational molecular lines have 
been detected in only four absor ption systems (with redshifts 0.25- 
0.89: IWiklind & Com bes' 1999, and references therein), despite in- 
tensive searches (e.g.fcurran et al. 2004). Based on the incidence of 
these lines and the number of continuum sources with fluxes larger 
than 10 mjy at 90 GHz we estimate that the number of molecular 
absorption systems observable with ALMA will be -^50, with only 
5-10 of these showing narrow lines - not enough to overcome the 
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uncertainties due to peculiar motions. Hence we liave decided not 
to pursue the case for molecular absorption lines any further. 



4 SENSITIVITY OF THE LYq FOREST TO RADIAL 
VELOCITY SHIFTS 

In Section[2]we have seen that the redshift drift is a very small ef- 
fect. In order to detect it, an experiment must achieve an overall ac- 
curacy with which radial velocity shifts can be determined of order 
~1 cm s~^ . In this section we will investigate how the properties of 
the Lya forest translate to a radial velocity accuracy, a^, and how 
depends on the instrumental characteristics of the spectra and on 
redshift. Specifically, we would like to know how many Lya forest 
spectra of which resolution and S/N are needed at a given redshift 
to achieve the required a-v ■ 

We will investigate these issues using artificial spectra. High- 
resolution observations have demonstrated that, to first approxima- 
tion, the Lya forest can b e decomposed in to a collection of individ- 
ual absorption lines (e.g. lKim et alJllOOll) . These are usually taken 
to be Voigt profiles and so each line is characterised by three pa- 
rameters: redshift, z, Hi column density, A'^hi, and velocity width, 
b. Here, we will reverse this decomposition process and generate 
(normalised) spectra with the desired instrumental characteristics 
from given lists of absorption lines. We will use two types of line 
lists. First, we will generate line lists from Monte Carlo (MC) sim- 
ulations based on the statistics of the largest available samples of 
absorption lines. Secondly, to validate our simulations, we will use 
8 line lists available in the literature that have previously been de- 
rived from high-resolution observations. 



4.1 Simulated absorption line lists 

We form simulated MC line lists by simply randomly drawing val- 
ues for the absorption line parameters fro m their observed distribu- 
tions (e.g. Huetal. 1995: JLu et alJl 19961 : Ixirkman & Tvtlejl 19971 : 
[Kim etal.,, 1997. , ,200 L ,20ol : 



Table 1. Observed hya forest line lists from the literature. 



exp 



2al 



(V) 



where 7 = 2.2, /3 = 1.5, 6 = 30 km s~^ and a;, = 8 km s'^ We 
impose limits of 15 < fe < 100 km s^^ and also restrict A^hi to 
the classical Lya forest regime, excluding Lyman limit and DLA 
systems: 12 < log A'^Hi(cm~^) < 16. The above distribution is 
normalised to give 10^ lin es with 13 64 < log AfHi(cm"^) < 16 
per unit redshift at 2; = 2 jKim et aLll200lh . The actual number of 
absorption lines in a given line list is drawn from a Poisson distri- 
bution with a mean determined by the normalisation. 

The above MC approach allows us to quickly generate large 
amounts of spectra with realistic characteristics. Note, however, 
that we make no assumptions regarding the underlying physics of 
the IGM in which the absorption occurs. We simply use the obser- 
vational fact that Lya forest spectra can be well represented as a 
random collection of Voigt profiles. 

The most significant difference between our MC line lists and 
the real Lya forest is clustering: the real Lya forest is not randomly 
distributed in redshift but shows significant r edshift-space correla- 
tions on scales of at least 100 km s~^ (e.g. ICristiani et alj 1 19951 : 
iFemandez-Soto et al.lll996l : lLiske et ai]|2000h . The impact of clus- 
tering will be discussed in detail in Section l431 



QSO 



2QSO 



A range" 



Reference 



QllOl-264 


2.145 


3226- 


-3810 


290 


1 


J2233-606 


2.238 


3400- 


-3850 


226 


2 


HE1122-1648 


2.400 


3500- 


-4091 


354 


1 


HE2217-2818 


2.413 


3550- 


-4050 


262 


2 


HE1347-2457 


2.617 


3760- 


-4335 


362 


1 


Q0302-003 


3.281 


4808- 


-5150 


223 


1 


Q0()55-269 


3.655 


4852- 


-5598 


535 


1 


Q0()00-26 


4.127 


5380- 


-6242 


431 


3 



Refere nces: 1 = lKim et alj i2002h . 2 = lKim et alj <200lh . 3 = ILu et alJ 

"Wavelength range covered by the Lya forest line lists in A. 
''Number of Lya forest lines. 



4.2 Real absorption line lists 

We have collected 8 QSO absorption line lists from the literature 
(see Table [T]!. These were derived from UVESA^LT (7 objects) 
or HIRES/Keck data (1 object). All spectra have a resolution of 
FWHM ~ 7 km s^^, while the typical S/N per pixel varies from 
~10 for QOOOO-26 to ~50 in the case of HE2217-2818. In all 
cases the absorption line lists were derived by fitting Voigt profiles 
to the spectra using VPFIT0 Details of the data acquisition and re- 
duction, as weU as of the line fitt ing a nd identification processes are 
given bvl Kim etalJ ( I2OOII . 12002*) and' Lu et al.l \\996). 

We note that two of the spectra, those of Ql 101— 264 and 
QOOOO— 26, contain DLAs, which effectively block out parts of the 
spectra. However, this is only expected to have a significant effect 
at high redshift, where the Lya forest line density is high, and so 
we have excluded the affected spectral region only in the case of 
QOOOO-26. 



4.3 The second epoch 

Simulating a i measurement requires a second epoch 'observation' 
of the same line of sight: we generate a second epoch absorption 
line list from the original line list (both real and simulated) by sim- 
ply shifting the redshift of each line according to a given cosmo- 
logical model: 



AZi = z{zi\ Ho, i^M, f^A) Ato, 



(8) 



where Ato is the assumed time interval between the first and sec- 
ond epochs. Unless stated otherwise, we use the standard, general 
relativistic cosmological model, assuming fiducial parameter val- 
ues of Ho = 70 hro km s"^ Mpc"\ = 0.3 and Oa = 0.7. 



4.4 Generating spectra 

Given an absorption line list a normalised spectrum is generated by 



(9) 



where A is the observed wavelength, A^ai is the number of absorp- 
tion lines in the spectrum, r is the optical depth of a Voigt pro- 
file and \a = 1215.67 A is the rest wavelength of the Hi Lya 



by R.F. Carswell et al, see |http://www.ast.cam. ac.uk/~rfc/vpfit.html] 
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Figure 5. (a) Example of an artificial Lya forest spectrum at 2 3 gen- 
erated from a Monte Carlo absorption line list, (b) Close-up of the region 
around 5000 A. (c) Artificial spectrum generated from the observed line list 
of Q0302-003. Both spectra have S/N = 100. 



transition. This spectrum is then pixelised, using a pixel size of 
0.0125 A, and convolved with a Gaussian line-spread function. Un- 
less stated otherwise, we will use a resolution element four times 
the pixel size, corresponding to a resolution of _R = 100 000 at 
5000 A. We then add random noise to the spectrum assuming Pois- 
son statistics, i.e. we assume a purely photon-noise limited exper- 
iment. All S/N values quoted in this paper refer to the S/N per 
pixel in the continuum. To begin with, we will only consider the 
spectral range between the assumed background QSO's Lya and 
Ly/3 emission lines, i.e. the classical Lya forest region. By con- 
struction this region cannot contain any H I transitions that are of 
higher order than Lya, and it covers an absorption redshift range of 
A/3/Aa(l + ZQso) — I < z < ZQso, where zqso is the QSO's 
redshift and Xp — 1025.72 A is the rest wavelength of the Hi 
Ly/3 transition. In Sections [5.1| and [J!2l we will also consider other 
spectral regions. 

In Fig.[5]we show examples of artificial spectra generated from 
both simulated (panels a and b) and real line lists (panel c). The ex- 
pected flux difference between two spectra of the same absorption 
lines, taken a decade apart, is shown in Fig. [6] A similar plot for 
the optical depth difference reveals that the variations discussed in 
Section [3.1.3l are unproblematic. 



4.5 Defining (t„ 

In principle, a z determination will involve the measurement of ra- 
dial velocity differences between the corresponding features of a 
pair of spectra of the same object taken several years apart. We 
now need a method to estimate the accuracy to which these differ- 
ences can be determined. This requires knowledge of how exactly 
the measurement will be performed. However, a priori it is not at 
all obvious what the optimal signal extraction method might be. To 
proceed nevertheless we choose to base our analysis on the generic 
concept of the total radial velo city information content of a sp ec- 
trum, which was developed bv iBouchv! Pepe & Oueio3 ( 1200 ih in 



"m. "a = 1-0. 0-0 
= 0.3, 0.0 

-n„. n, = 0.3, 0.7 




= 22 cm s 
Av = 8.4 cm s~ 
Av = 2.8 cm 
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Figure 6. Flux difference between two ailificial, noiseless spectra of the 
same Lyo forest at 2 3 simulated for two observing epochs separated by 
Ato = 10 yr and for various combinations of and Q/^ as indicated. 
The redshift drift implied by these pai'ameters is also given. 



the context of optimising radial velocity searches for extra-solar 

planets. 

Following Bouchv et al. I jlOOlh we begin by expressing the 
flux observed in pixel i at the second epoch as a small perturbation 
on the first epoch flux in the same pixel: 

S2r =Su + ^— A., (10) 

dA c 

which defines a small velocity shift Avi for each pixel. \i is the ob- 
served wavelength of the ith pixel and dSi /dA is the spectral slope 
of the flux at that pixel. To first order the slope does not change 
between the two epochs and hence it carries no epoch designation. 
Averaging the velocity shift over all pixels in a spectrum, using 
weights Wi, we have: 



Aw : 



(11) 



Clearly, the weight for the ith pixel should be chosen as the inverse 
variance of Avj. In calculating this variance we must differ from 
iBouchv et alj feOOlh . In the case of stars one of the spectra can 
be assumed to be a perfect, noiseless template, essentially because 
additional information on the same type of star can be used to define 
it. However, since every Lya forest spectrum is unique we cannot 
make the same assumption here, so that in our case both spectra 
have noise. Hence we find: 



A, 



dSi 



2 



2 

+ (72i 



+ 



(S2i ~ Slif' 2 
V dA j 



(12) 



where a\i and a2i are the flux errors in the ith pixel of the first and 
second epoch spectra, respectively, and a^i is the error on the slope 
of the flux at pixel i. We can see that a low weight is assigned to 
noisy pixels and to those that have a small gradient, i.e. pixels in the 
continuum or in the troughs of saturated absorption lines. Finally, 
with the above choice of weights the error on Av is given by: 



(13) 



The above process has the advantage of conveniently attach- 
ing a single figure of merit to a given pair of spectra in a non- 
parametric, model-independent way: (7„ simply represents the fun- 
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damental photon-noise limit of the accuracy to which an overall 
velocity shift between the two spectra can be determined. It is es- 
sentially just a measure of the 'wiggliness' and of the S/N of the 
spectra. However, we point out that cr„ does not entirely capture 
all of the information contained in a pair of spectra with respect 
to a i measurement. From Fig. [2] it is clear that the difference be- 
tween the first and second epoch spectra is not simply an over- 
all velocity shift. The second epoch spectrum will also be slightly 
compressed with respect to the first epoch spectrum because the 
redshift drift is larger at higher redshifts (i.e. longer wavelengths) 
than at lower redshifts. The additional information that is contained 
in this small alteration of the spectrum's shape is entirely ignored 
by the above method [because of the simple averaging operation in 
equation dllH and hence it is clearly a sub-optimal method of es- 
timating the sensitivity of a pair of spectra to z. Despite this short- 
coming we will retain the above definition of CTi, for the sake of its 
simplicity. 



4.6 Results 

We are now ready to derive the relevant scaling relations for (7„ . We 
begin by making two points regarding equation il2\ . First, we note 
that for a fixed total integration time (the sum of the integration 
times spent observing the first and second epoch spectra) the sum 
+ in equation l ll2t takes on its minimum when the first 
and second epoch integration times are equal. Hence, the smallest 
possible cr„ is only achieved when the spectra of both epochs have 
the same S/N. In the following we will assume that this is the case. 
Secondly, from equation il2i it is clear that av scales as (S/N)^^, 
as expected for a photon-noise limited experiment. 

Consider now a set of Nqso targets that all lie at the same 
redshift, each of which has been observed at two epochs such that 
all 2A*'qso spectra have the same S/N. Again, since we are consid- 
ering a purely photon-noise limited experiment, ct„ should scale as 

— 1/2 

-^pix (where A^pix is the total number of independent data points 

— 1/2 

in the sample) and hence also as Nqqq . Furthermore, in this ideal 
case, it is irrelevant how the total S/N is divided among the targets. 
Hence, for simplicity we will use A'^qso = 1 in the following. 

We now examine the behaviour of (t„ as a function of redshift 
using the MC absorption line lists. For various QSO redshifts in 
the range 2 ^ zqso ^ 5 we have generated 10 pairs of Lya forest 
line lists and spectra as described in Sections [4.1| to |4.4l where each 
spectrum was given a S/N of 13 000. We then measured each pair's 
cr„ according to equation J13b . The result is shown as blue dots in 
Fig. |7] where each point and error bar is the mean and ±1 r.m.s. 
of the 10 individual measurements at each redshift. We stress that 
we are plotting the expected accuracy of a velocity shift measure- 
ment performed on a single pair of spectra of a single target at a 
given redshift, where each spectrum has S/N = 13 000. We are not 
plotting the combined accuracy of 10 such pairs. 

From Fig.Qwe can see that the radial velocity sensitivity im- 
proves rather rapidly with redshift for zqso < 4, but the decrease 
is somewhat shallower at zqso > 4. Overall, (t„ improves by a 
factor of almost 3 when moving from zqso = 2 to 5. Specifically, 
we find: 



(1 
(1 



ZQSOj 

zqso) 



-1.7 
-0.9 



ZQSO < 4 
ZQSO > 4 



(14) 



• MC line lists 
A 7 = 

X 7 = 0, Az = 0.4 




The behaviour of (t„ as a function of redshift is due to the combina- 
tion of several factors. The first is the redshift evolution of the Lya 
forest line density (cf. equation O. At higher redshift more spec- 



Figure 7. The blue dots with en'or bars show the accuracy with which a ra- 
dial velocity shift can be determined from a hya forest spectrum as a func- 
tion of QSO redshift. The red line pai'ameterises the redshift dependence as 
in equation 114) . Each point is the mean cr^ measured from 10 pairs of arti- 
ficial spectra with S/N = 13 000, generated from simulated MC absorption 
line lists. The error bars show the ±1 rm.s. range of the 10 simulations. The 
green triangles show the results for simulations where the redshift evolution 
of the hya forest has been switched off, i.e. where 7 = (cf. equation|7). 
The orange crosses show the result of also restricting the cr„ measurement 
in each spectrum to a redshift path of constant length A 2 = 0.4, as op- 
posed to using the full hya forest region between the QSO's hya and Ly/3 
emission lines. 



tral features are available for determining a velocity shift and so 
cr„ decreases. However, from z ~ 4 the absorption lines severely 
blanket each other and the number of sharp spectral features does 
not increase as rapidly anymore, causing the flattening of CTj, at 
Zqso > 4. The green triangles in Fig. |7] show the a-v measure- 
ments that result from simulations where the redshift evolution of 
the Lya forest has been switched off, i.e. where the evolutionary 
index 7 has been set to (cf. equation |7). Indeed, we can see that 
in this case there is no evidence of a break. 

Secondly, we recall that each Lya forest spectrum covers the 
entire region between the QSO's Lya and Ly/3 emission lines. The 
redshift path length of this region is given by Az — 0.156(1 + 
2qso)- Hence the number of independent pixels per spectrum also 
increases as (1 -I- zqso)- Since the S/N per pixel is kept constant 
this implies a larger number of photons per spectrum and hence an 
improved sensitivity to radial velocity shifts. The effect of this can 
be seen by comparing the green triangles in Fig.|7]with the orange 
crosses which are the result of using a constant redshift path length 
of Az = 0.4 for each cr„ measurement, as well as 7 = 0. 

Finally, with 7 = and Az = const, the sensitivity to wave- 
length shifts should be constant as a function of zqso, and so the 
sensitivity to velocity shifts should go as (1 + zqso)^^ (as can be 
seen from equation II 2t. In fact, the orange crosses in Fig. [7] de- 
crease more slowly than this because the widths of the absorption 
lines in wavelength space increase as (1 + zqso), making the edges 
of the lines less steep and hence slightly decreasing the spectrum's 
sensitivity to wavelength shifts. 

In Fig. [8] we compare these results derived from the MC line 
lists to those from the real line lists. From Table [T] we can see that 
the real line lists (and hence the corresponding spectra) do not cover 
the full LyQ forest regions, with differently sized pieces missing 
both at the low and high redshift ends. Therefore we must correct 
the Ov values derived from the real line lists in order to make them 
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Figure 8. Comparison of ctv measurements derived from simulated MC line 
lists (blue dots and solid line, same as in Fig.|2) and from real line lists (red 
stars). The green circles show the results derived from simulated MC line 
lists that include a simple scheme for clustering absorption lines in redshift 
space (see text for details). 



directly comparable to the values from the simulated lists. The cor- 
rection is achieved by first assigning a new, slightly different QSO 
redshift to each spectrum, such that the 'missing' low and high red- 
shift parts of the Lya forest region are equally large. We then de- 
crease the measured a„ by a factor (Azobs/Az)^''^, where A^obs 
is the redshift path length covered by the observed line list and Az 
is the redshift path length of the full Lya forest region at the new 
zqso- The correction factors range from 0.56 to 0.99. 

The red stars in Fig. [8] show the corrected values derived 
from single pairs of spectra generated from the real absorption line 
lists with S/N = 13 000, while the blue dots show the measure- 
ments from the MC line lists (same as in Fig.[7](. Overall the agree- 
ment between the results from the MC and real line lists is very 
good, particularly at high redshift. At zqso ~ 2.4 the (t„ values 
from the real line lists are ~15 per cent higher than those from the 
MC lists. By far the most significant deviation occurs at the lowest 
redshift where the of QllOl— 264 is higher than expected by 
47 per cent. However, this is not too surprising as the line of sight 
towards QllOl— 264 is known to pass through an unusually low 
number of absorbers with log Ai'Hi(cm~^) > 14 jKim et aUl2002h . 

We believe that the small differences at zqso ~ 2.4 are 
mainly due to clustering of real absorption lines in redshift space. 
Clustering has the effect of reducing the number of spectral features 
because it increases line blanketing. However, at high redshift line 
blanketing is already severe because of the high line density and 
so clustering has a relatively smaller effect at high redshift than at 
low redshift. We demonstrate that clustering can explain the ob- 
served differences by generating a new set of MC line lists which 
incorporate a toy clustering scheme: first, we randomly draw the 
positions of 'cluster' centres from the Lya forest redshift distribu- 
tion. We then populate each 'cluster' with n ^ absorbers, where 
n is drawn from a Borel distribution l lSaslawll"9^8ft ). Since n can be 
0, 1 or > 1 this process generates voids, single 'field' absorbers, 
as well as groups and clusters of lines. Finally, absorption lines are 
distributed around their 'host' clusters according to a Gaussian dis- 
tribution with FWHM = 120 km s^^. The cr„ values that result 
from this new set of MC line lists are shown as open green circles 
in Fig. [8] The increase of Uu compared to the unclustered simu- 
lations is clearly very similar to that observed for the real line lists 
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Figure 9. Radial velocity accuracy as a function of spectral resolution. Each 
point and error bar is the mean and ±1 r.m.s. of 10 simulations at zqgo = 
4 and S/N = 13 000. The pixel size is kept constant and is chosen such 
that a resolution element is sampled by 3 pixels at the highest resolution 
considered. The labels along the bottom axis denote the resolving power at 
5000 A, while the labels along the top axis show the equivalent FWHM of 
an unresolved line. The vertical red line marks the mode of the absolution 
lines' b parameter distribution (cf. equation|7). 



and we conclude that clustering can indeed explain the small differ- 
ence between the results obtained from the observed and simulated 
line lists at zqso ~ 2.4. In any case, the near coincidence of the 
cjv value of J2233— 606 at zqso = 2.23 with the expected value 
demonstrates that not all lines of sight are adversely affected by 
clustering. In the following we will assume that such sightlines can 
be pre-selected and hence we will ignore the effects of clustering 
in the rest of this section. 

We now turn to the behaviour of (j„ as a function of spectral 
resolution. For various resolving powers in the range 2500 i? 
1.33 X 10^ we have generated 10 pairs of line lists and spectra with 
zqso = 4 and S/N = 13 000, and measured their values as be- 
fore. The result is presented in Fig.|9] where we show the resolving 
power along the bottom axis and the corresponding FWHM of an 
unresolved line along the top axis. We stress that the pixel size was 
the same for all spectra (— 0.0125 A) and that it was chosen such 
that a resolution element was well sampled even at the highest reso- 
lution considered. Hence the strong increase of a„ towards lower R 
in Fig.|9]is not due to different numbers of photons per spectrum or 
sampling issues. Instead, it is simply due to the loss of information 
caused by convolving the spectrum with a line-spread function that 
is broader than the typical intrinsic absorption line width (marked 
by the vertical line in Fig.|9). Indeed, at _R > 30 000 the Lya forest 
is fully resolved and in this regime is independent of R. 

Summarising all of the above we find that the accuracy with 
which a radial velocity shift can be determined from the Lya forest 
scales as: 



= 2 



/_S^\"^ / jVQSO \~^ / 1 + ZQSO \ - 

V2370y V 30 y V 5 / 



.(15) 



where the last exponent changes to —0.9 at zqso > 4, and where 
the same S/N (per pixel) is assumed for all Nqso spectra at both 
epochs. 
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5 SENSITIVITY GAINS FROM OTHER SPECTRAL 
REGIONS 

In the previous section we have investigated the sensitivity of the 
Lya forest to radial velocity shifts, using only the H I Lya transi- 
tion in the region between a QSO's Lya and Ly/3 emission lines. 
However, modern echelle spectrographs are capable of covering a 
much wider spectral range in a single exposure and so the question 
arises whether other spectral regions, containing absorption lines 
from other ions or other H I transitions, can expediently contribute 
towards a z measurement. 

Before applying the procedure of the previous section to more 
extended spectra, we briefly follow-up on the discussion at the end 
of Section [43] where we pointed out that does not capture the 
information that is contained in the i-induced change of the shape 
of a spectrum. In the case of a Lya-only spectrum this change of 
shape consisted of a compression of the spectrum. If we allow addi- 
tional transitions with different rest wavelengths then we no longer 
have a one-to-one correspondence between absorption redshift and 
wavelength, and so the redshift drift will in general induce a much 
more complex change of the shape of a spectrum. Hence it is clear 
that any attempt to harness this additional information must involve 
the complete identification and modelling of all absorption features 
used in the analysis. 

Indeed, the complete identification of all metal absorption 
lines will be necessary in any case, even if one endeavours to mea- 
sure the redshift drift only from the Lya forest. The point is that the 
Lya forest may of course be 'contaminated' by metal lines. Since 
these lines may arise in absorption systems that lie at completely 
different redshifts compared to the Lya lines, their redshift drift 
may also be very different. Hence, any lines wrongly assumed to 
be Lya could potentially lead to erroneous results. The only way to 
avoid such biases is to completely identify all absorption features. 



5.1 The Ly/3 forest 

Each H I absorption line in the Lya forest has corresponding coun- 
terparts at shorter wavelengths that result from the higher order 
transitions of the Lyman series. Obviously, all the arguments con- 
cerning the suitability of the Lya forest for a i measurement also 
apply to these higher order lines, and so, in principle, they should 
be almost as useful as the Lya lines (where the qualifier 'almost' is 
owed to the decreasing optical depth with increasing order). How- 
ever, in practice one has to contend with (i) confusion due to the 
overlap of low-redshift, low-order lines with high-redshift, high- 
order lines, and (ii) with an increased uncertainty in the placement 
of the QSO's continuum. Both of these difficulties are aggravated as 
one proceeds up the series towards shorter wavelengths. For these 
reasons we will not consider any H I transitions that are of higher 
order than Ly/3. 

We now extend our simulated spectra by adding the region be- 
tween a QSO's Ly/3 and Ly7 emission lines immediately bluewards 
of the Lya forest. This region extends the redshift path for Lya 
lines by a factor of 1.28 towards lower redshifts and also contains 

Ly /3 lines in the redshift range At/ A^(l + 2QSo) — 1 < z < zqso, 
where A^. = 972.54 A. Note, however, that the real Lya line lists 
do not extend below the QSOs' Ly/3 emission lines (cf. Table [Til. 
Hence we are forced to supplement the real Ly/3 absorption lines in 
this region with Lya lines drawn from the MC simulations. 

In Fig.[lO]we compare the (j„ measurements derived from the 
extended MC line lists (blue dots) to the Lya-only results (solid 
line, same as in Figs.|7]and[8j. We find that including the Ly/3 for- 
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— MC line lists: Lya only 
• MC line lists; Lya + Ly^ 
real line lists: Lya only 
real line lists: Lya + Ly|S 
□ real line lists: Lya + LyjS + metals 
A reol line lists: full spectrum 




Figure 10. Comparison of av measurements derived from Lyo-only line 
lists and various extended line lists covering additional spectral regions and 
including other absorption lines. The solid line and solid stars show, respec- 
tively, the results from the Lyo-only MC and real line lists as in Fig. [8] The 
blue dots and open stars show the corresponding improved (t„ values that 
result from the addition of the Ly/3 forest. For the real line lists, the open 
squares show the effect of further adding the available metal lines in the 
Lycn-Ly/3 forest region. Finally, for two of these QSOs, we show as solid 
triangles the outcome of using all lines accessible in existing spectra. 



est improves cr„ by a factor of 0.88 ± 0.006. The corresponding 
comparison for the real line lists (open and solid stars) yields im- 
provement factors of 0.82 to 0.92, with a mean of 0.88, in very 
good agreement with the MC results. 



5.2 Metal lines 

In addition to the H I absorption, every high-redshift QSO spectrum 
shows absorption lines from a number of other ions, such as C IV, 
Si IV or Mg II. Although these metal lines are far less numerous 
than the H I lines, they are also much narrower: their widths are of 
order a few km s^ ^ . In fact, many metal lines are unresolved in cur- 
rent spectra and so the widths of the narrowest lines are unknown. 
This suggests that metal lines may considerably increase a spec- 
trum's sensitivity to radial velocity shifts and may hence supply 
valuable additional constraints on i. 

Could peculiar motions vitiate this supposition? This ques- 
tion is difficult to answer with certainty because the exact ori- 
gin of many of the various classes of metal absorption lines is 
still under debate. It is nevertheless evident that the structures 
responsible for the absorp tion form an inhomogeneous set (e.g. 
IChurchill et al.l2000l[2007l) and represent a large range of environ- 
ments, depending on the absorbing ion, column density and red- 
shift. For example, low column density CiV absorbers probe the 
IGM (e.g. ISimcoe et al.ll2004l : [Songai la"2Q06') while strong Mgll 
lines are associated with galaxies ( e.g. Bergeron & Boisse 199l|; 
ISteidel. Dickinson & Perssonlll994l : Ebetti et alj|2005i) . However, 
we note that absorbers associated with galaxies are generally found 
to have absorption cross-sections on the order of tens of kpc (e.g. 
ISteidellll995l : IChurchill et al] l2OO0l: lAdelberger etaLl |2005|). For 
such distances and for halo masses of ~10^^Mq dSteidel et al.l 
1 19941 ; iBouche et alj |2004 l2006l) one derives accelerations of a 

few X 10~* cm s~^, while absorber kinematics point to pe- 

-1 



culiar velocities of no more than a few hundred km s ^ (e.^ 
IChurchill et"ai]| 19961) . According to Appendix [A] peculiar motions 
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Table 2. Observed lists of metal absorbers. 



Name A range" Reference 



Ql 101 -264 


3050- 


-5765, 


225 


this work 




5834- 


-8530 






J2233-606 


3400- 


-3850 


49 


2 


HEl 122- 1648 


3500- 


-4091 


18 


1 


HE2217-2818 


3550- 


-4050 


59 


2 


HE1347-2457 


3760- 


-4335 


48 


1 


Q0302-003 


4808- 


-5150 


5 


1 


Q0055-269 


4852- 


-5598 


14 


I 


QOOOO-26 


4300- 


-6450, 


100 


this work 




7065- 


-8590, 








8820- 


-9300 







References: 1 = lKim et alj j2002h . 2 = lKim et alj i200ll) . 

"Wavelength range(s) searched for metal lines in A. 
''Number of metal absorbers, defined as the number of 
unique sets of {ion, z, N, b}. 



of this magnitude are not problematic. The main caveat here is the 
unknown fraction of strong metal absorbers arising in starburst- 
driven outflows (e .g. lAdelbergerer^l2003l : Isimcoe et al.ll2006l : 
iBouche et al.l200^ . see also Section [3.1.2t which must be expected 
to experience much larger accelerations. Howe ver, it may be possi- 
ble to identify these systems kinematically (P rochaska et al ]|2008l) 
and hence to exclude them from the analysis. In any case, in this 
section we will proceed on the assumption that peculiar motions do 
not generally invalidate the use of metal lines for a i measurement. 

Equation Q gave a succinct parameterisation of the proper- 
ties of the Lya forest. Unfortunately, equivalents for all the various 
metal line species do not exist, and hence we are unable to gauge 
the effects of metal lines usi ng M C simula t ions. 

iKim et alj ( l200lLl2002h and lLuetalJ ^199^. from whom we 
obtained the real Lya forest line lists used in the previous sections, 
also identified metal lines in their spectra. Unfortunately, they only 
published measured parameters for lines lying in the Lya forest 
region. Although these lines will also have transitions elsewhere 
(which were probably used in the identification process) the pub- 
lished line lists do not provide a complete view of the metal line 
population outside of this region. 

To improve on this situation we went back to the spectra of two 
of our QSOs, Ql 101-264 and QOOOO-26, and derived our own 
metal line lists covering the entire spectral range available to us. In 
the case of Q 1 1 1 — 264 the spectrum was the same as that used by 
iKim et al.1 ilOOA whereas for QOOOO-26 we used a UVES/VLT 
spectrum of similar quality as the HIRES/Keck spectrum used by 
iLu et al. |(l996). In the following we will refer to these two line lists 
as the 'complete' metal line lists. 

Table|2]summarises the wavelength ranges that were searched 
for metal lines and the total number of metal absorbers that were 
found in each QSO spectrum. The simulated spectra are generated 
in exactly the same way as in Section lSTT] except that we now add 
all available metal lines, and extend the spectra of Ql 101 — 264 and 
QOOOO-26 to the red limits given in Tabled 

Fig. 1 101 shows the effect of the metal lines on a^- The green 
open squares show the result of only using the Lya-l-Ly/3 forest re- 
gion as in Section [5TT] but adding in all the available metal lines in 
this region. Comparing this to the H 1-only results of the previous 
section (open red stars) we find that the metal lines improve 
by factors of 0.58 to 0.99, with a mean of 0.85. We point out that, 
strictly speaking, the derived (j„ values are only upper limits for the 



six line lists taken from the literature because of their incomplete 
coverage of the Ly/3 forest region. Indeed, the best improvement is 
achieved for one of the complete line lists, Ql 101— 264 (the low- 
est redshift QSO), which has a particularly rich metal absorption 
spectrum (cf. Table|2]l- 

Finally, for our two complete line lists we show as solid green 
triangles the effect of also adding in the accessible spectral regions 
redwards of the Lya forest. The additional metal lines further im- 
prove the a„ values of Q 1101-264 and QOOOO-26 by factors of 
0.55 and 0.90, respectively. With only two values it is obviously 
impossible to draw a firm conclusion regarding the average im- 
provement offered by the metal lines redwards of the Lya forest. 
We therefore choose to be conservative and adopt the larger of the 
two as a typical value. 

Summarising Sections 15.11 and 15.21 the above experiments 
have shown that the normalisation of equation l |15l l can be reduced 
by a factor of 0.88 x 0.85 x 0.9 = 0.67 by considering not just 
the Lya forest but all available absorption lines, including metal 
lines, over the entire accessible optical wavelength range down to a 
QSO's Ly7 emission line. Hence we now obtain: 



6 MULTIPLE EPOCHS 

Fundamentally, a redshift drift experiment consists of simply mea- 
suring the velocity shift between two spectra of the same QSO(s) 
taken at two distinct epochs separated by some time interval Ato- 
This is the view we took in the previous two sections where we de- 
termined the fundamental photon-noise limit of the accuracy with 
which this shift can be measured, and its scaling behaviour. How- 
ever, in practice this notion is too simplistic. First of all, it implicitly 
assumes that the total integration time, tint, required to achieve the 
necessary S/N would be negligible compared to Ato, so that the 
two epochs of observation are well-defined. As we will see in the 
next section this assumption is not valid. Secondly, for a variety of 
reasons it may be desirable to spread the observations more evenly 
over the whole period Ato instead of concentrating them in just two 
more or less well-defined epochs at the interval's endpoints. Hence 
the question arises how the accuracy of a redshift drift experiment 
is affected by distributing the total available tint over multiple ob- 
serving epochs within the interval Ato- 

Let us assume then that observations take place at Ne different 
epochs, where the jth epoch is separated from the first by Atj, so 
that Ati = and AtM^ ~ Ato. We can straightforwardly gener- 
alise the framework developed in Section |43] bv turning equation 
dlOl ) into a continuous equation for the expected normalised flux in 
the ith pixel at time At: 

S^iAt) = So,: + ^ At = So, + m,At. (17) 

dA c 

The idea is now to fit this linear relation to the observed fluxes 
Sji at times Atj with errors Oji, yielding an estimate of the slope 
rrii and hence of iii for each pixel. (Note that Soi is a nuisance 
parameter. It represents the 'true' flux at the first epoch as opposed 
to the observed value Su.) The maximum likelihood estimator for 
rrii is 



At2 - At 

where the bar denotes the weighted average over all epochs: 
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The variance of irii is given by 



J2 "^-^ ( 



Af2 - /\t' 



(19) 



(20) 



With rui and its variance in place we can write down the equivalent 
of equation (12}: 



+ 



(21) 



which in turn allows us to compute ii averaged over all pixels and 
its error, ai, corresponding to equations dl lb and ( I13t . Finally, we 
re-define (j„ = Ato. This new version of cr^ now includes the 
effect of multiple observing epochs and an arbitrary distribution of 
the total integration time among them. It is straightforward to show 
that for Ne — 2 one recovers exactly the original av of Section l431 
We are now in a position to amend equations dlSt and il6\ 
to include the effect of multiple epochs. Recall that these scaling 
relations were derived for the case of A^'e = 2 and for equal S/N 
in the spectra of both epochs. Since the variance of the normalised 
flux in pixel i scales as the inverse of the integration time we now 



write an as 



as 



(22) 



where fj is the fraction of the total tint used at the jth epoch 
fj ~ 1), and (Ji denotes the flux error (in the ith pixel) that 
one would obtain if half of the total tint were used. Further defining 



Atj = hjAto, we can re- write equation i20\ as 



= al^{N, = 2,/i = /2 = 0.5) /(Are,/ii...iVe,,fi...ivJ.(23) 

The first term above is just the variance of rrii that one obtains in 
the case of A^'e = 2 and equal splitting of tint. The second term is 
a 'form factor' that only depends on the distribution of tint within 
Ato. Again, it is straightforward to show that g{Ne = 2, /i — 
f2 = 0.5) — 1. Since the form factor is the same for all pixels, and 
since ct^. is the dominant term in equation l l21t . the sought-after 
modification of the scaling relation amounts to simply applying 
the form factor: 



= 1.35 
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(24) 



Note that the symbol 'S/N' now refers to the total S/N per object 
accumulated over all epochs, in contrast to equations dlSt and ( I16t 
where it referred to the S/N achieved in each of two epochs. 

Note also that we have dropped the dependence of the form 
factor on hi...N^ by considering every night within the period Ato 
as a potential epoch. This fixes Ne and hj = (j — l)/{Ne — 1), 
while fj is constrained to lie in the range ^ fj ^ l/ti^t, where / 
is the length of a night (which we will assume to be 9 h on average). 
Thus we find: 



<7(iVe,/l...ivJ = 




.(25) 



Clearly, fj will be for most nights. Nevertheless, there are 
obviously a large number of different possible distributions for 
the fjS. The best distributions are those that are symmetric and 
peaked towards the endpoints of Ato. A flat distribution, with equal 
observations taking place on n equally spaced nights, results in 
g = y/3 (n- l)/(n + 1) ^ 1.7 for n ^ 1. Thus the otherwise 
quite desirable arrangement of observing at a more or less constant 
rate throughout the period Ato comes with a rather severe penalty 
attached. A priori, it is difficult to estimate the best g value that 
can be realistically achieved in practice. From now on we will as- 
sume, perhaps somewhat arbitrarily, that all observations occur as 
much as possible towards the beginning and end of Ato with the 
constraint that the observing rate averaged over some intermediate 
timescale of, say, a month cannot exceed 1/3, i.e. that no more than 
a third of any month's telescope time is used for the redshift drift 
experiment. Depending on the ratio of tint and Ato this results in g 
values of ~1.1. Essentially, this configuration simply shortens the 
effective length of the experiment by the amount of time it takes to 
complete the observations at either end of Ato. 



7 CAN WE COLLECT ENOUGH PHOTONS? 

In Sections|4l[5]and[6]we learnt what S/N ratio is required to achieve 
a given sensitivity to radial velocity shifts using QSO absorption 
spectra. In a photon-noise limited experiment the attainable S/N de- 
pends only on four quantities: the brightness of the source, the size 
of the telescope's collecting area, the total integration time and the 
total efficiency. By 'total efficiency' we mean the ratio of the num- 
ber of detected photo-electrons to the number of source photons at 
the top of the atmosphere, i.e. it comprises atmospheric absorption 
and all losses occurring in the combined telescope/instrument sys- 
tem, including entrance aperture losses and the detector's quantum 
efficiency. 

In this Section we will investigate in detail the 5-dimensional 
parameter space that is spanned by the above four quantities and 
redshift, in order to determine whether a feasible combination ex- 
ists that would allow a meaningful z measurement. 



7.1 S/N formula 

We begin by writing down the relation between the S/N per pixel 
and the above four parameters for the photon-noise limited case: 



S^ 
N 



700 



Zx_ 
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10 h 0.25 



(26) 



where D, tint and e are the telescope diameter, total integration 
time and total efficiency, Zx and mx are the zeropoint and appar- 
ent magnitude of the source in the X-band, respectively, and Zr — 
(8.88 X 10^°) s-i m"^ nm'^ is the AB zeropoint ( IOkeiiT9 74) for 
an effective wavelength of 6170 A [corresponding to the Sloan Dig- 
ital Sky Survey (SDSS) r-band]. The normalisation of the above 
equation assumes a pixel size of 0.0125 A (see Section|4j4) and a 
central obscuration of the telescope's primary collecting area of 10 
per cent. D — 42 m correspon ds to the Baseline Reference D esign 
for the European ELT (E-ELT: lGilmozzi & Spvromilioll2007h . 



7.2 High-redshift QSOs 

The photon flux from QSOs is of course not a free parameter that 
can be varied at will. Instead we will have to content ourselves 
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with what will be offered by the population of real QSOs known 
at the time of a hypothetical i. experiment. Here we do not wish to 
speculate on possible future discoveries of QSOs and hence we will 
restrict ourselves to the ones known already today. In the following 
we will extract a list of potential targets for a i experiment from 
existing QSO catalogues. For each candidate target QSO we will 
need a reliable magnitude that can be used to estimate its photon 
flux, as well as its redshift. 

The largest QSO catalogue with reliable, homogeneous pho- 
tometry and redshifts curr ently available is the fourth edition of 
the SDSS Quasar Catalog l lSchneider et al.ll2007 t). Being based on 
the fifth data release of the SDSS, it yields 16 913 QSOs with 
zqso ^ 2. The catalogue provides PSF magnitudes in the ugriz 
bands which we do not correct for Galactic extinction (as is ap- 
propriate for S/N calculations). Since we are interested in the con- 
tinuum flux we will use, for each QSO, the magnitude of the 
bluest filter that still lies entirely redwards of the QSO's Lya emis- 
sion line. Specifically, for objects with zqso < 2.2 we will use 
the p-band magnitude; for 2.2 ^ zqso < 3.47 the r-band; for 
3.47 ^ ZQSO < 4.61 the i-band; and for 4.61 ^ zqso the z-band. 
We then apply a small correction to the selected magnitude to trans- 
form the observed flux to that expected at the centre of the Lya for- 
est assuming a power-law spectral shape of the form oc 
( lFrancislll993h . 

Unfortunately, the SDSS catalogue does not cover the whole 
sky. The largest collection of QSOs covering the entire sky is the 
12th edition o f the Catalogue of Quasars and Active Nuclei recently 
published bv I Veron^ Cettv & Veroiil j2006h , which contains many 
additional bright QSOs not included in the SDSS catalogue. How- 
ever, since the Veron catalogue is a compilation of data from many 
different sources its photometry is very inhomogeneous and cannot 
readily be converted to a photon flux. We will dispense with this 
inconvenience by matching the Veron catalogue to the all-sky Su- 
perCOSMOS Sky Survey (SSS: lHamblv et al.ll200ll) . Although the 
photographic SSS photometry is not endowed with reliable abso- 
lute calibrations either, at least it is homogeneous and covers three 
bands (_Bj, R and /), allowing us to synthesise approximate SDSS 
magnitudes. 

We proceed by first selecting all QSOs from the Veron cata- 
logue with ZQso J? 2 and not flagged as unreliable, resulting in 
21 895 objects. For each of these we then identify the nearest SSS 
object, allowing a maximum separation of 5 arcsec, and retrieve the 
corresponding SSS catalogue data. 938 Veron objects have no SSS 
counterpart, presumably because their coordinates are unreliable. 
We then use those 11 744 objects that are in common to the SDSS 
and combined Veron-SSS catalogues to derive linear relations be- 
tween the SDSS and SSS magnitudes, allowing for zeropoint off- 
sets and colour terms. Such relations are reasonable representations 
of the data and we find that the distributions of the residuals have 
r.m.s. values of less than 0.3 mag in all cases. Finally, we purge the 
common objects from the Veron-SSS catalogue and use the above 
relations to synthesise an SDSS magnitude in the appropriate band 
(see above) for each of the remaining QSOs. 

For those QSOs in the initial Veron catalogue that have no 
match in the SSS, or which are missing an SSS band needed to 
synthesise the required SDSS magnitude, we will simply use the V 
or i?-band magnitude as listed in the Veron catalogue, provided it 
is non-photographic. 

In summary, the final combined sample of 25 974 QSOs is 
constructed from three sub-sets: (i) SDSS; (ii) objects with redshifts 
from the Veron catalogue and photometry from the SSS (converted 
to the SDSS system); and (iii) objects where both the redshifts and 




Figure 11. The dots show the known, bright, high-redshift QSO population 
(separated by sub-sets as indicated, see text) as a function of redshift and 
estimated photon flux at the centre of the Lya forest. Along the right-hand 
vertical axis we have converted the photon flux to a corresponding Johnson 
V-band magnitude. The background colour image and solid contours show 
the value of CTi, that can be achieved for a given photon flux and redshift, 
assuming D = 42 m, e = 0.25 and ti^t = 2000 h. The contour levels 
are at = 2, 3, 4, 6, 8 and 10 cm s~^. The dotted contours show the 
same as the solid ones, but for D = 35 m or, equivalently, for e = 0.17 or 
tint = 1389 h. 



the photometry are taken from the Veron catalogue. We remind the 
reader that the quality and reliability of the photometry decreases 
rapidly from (i) to (iii). 



7.3 Achievable radial velocity accuracy 

In Fig.[TT|we plot our QSO sample, split by the above sub-sets, in 
the A'^phot-ZQso plane, where A^phot is a QSO's photon flux at the 
top of the atmosphere and at the centre of the QSO's Lya forest, 
as implied by the appropriate magnitude described above. Using 
equations l l24t and i26l and assuming values for D, e and tint we 
can calculate, for any given combination of A^'pi^ot and zqso^ the 
value of CTi, that would be achieved if all of the time tint were in- 
vested into observing a single QSO with the given values of A'^piiot 
and zqso- The background colour image and solid contours in 
Fig.ll llshow the result of this calculation, where we have assumed 
D = 42 m, e = 0.25 and tint = 2000 h. Note that we have in- 
cluded both the improvement of cjv afforded by the Ly/3 forest and 
the metal lines as well as the deterioration caused by spreading tint 
over a 0.9 yr period at either end of a Ato = 20 yr interval. 

From Fig. [TT] we can see that, although challenging, a rea- 
sonable measurement of z[z) is within reach of a 42-m telescope. 
There exist a number of QSOs that are bright enough and/or lie 
at a high enough redshift to provide reasonable values of (t„. We 
find 18 objects at a„ < 4 cm s~^ and 5 objects at < 3 cm s^^, 
with good coverage of the redshift range 2-4. One object even gives 
(Ji, = 1.8 cm s^^. However, for a smaller telescope with D = 35 m 
the number of objects with cr^ < 4 cm s~^ reduces to only 7 (cf. 
dotted contours). 

In Fig. [12] we show more comprehensively how the number of 
QSOs with Ov smaller than a given value depends on the telescope 
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Figure 12. The colour image and the contours show the number of QSOs 
for which the value on the ordinate or better can be achieved for a given 
combination of telescope size, efficiency and integration time. The contour 
levels are at Nqsq = 3, 5, 10 and 30. 

parameters and integration time, which we summarise into a single 
'normalised observational setup parameter', O: 

0=(^Y (27) 
V42m/ 0.25 2000h ^ ' 

(cf. equation[26ll. For a given value of O the colour image and con- 
tours in Fig. [12] show the number of QSOs that would give a (7„ 
equal to or smaller than the value along the ordinate if all of tint 
was spent on any one of them. For example, if we wanted to be 
able to choose our targets from 30 QSOs bright enough and/or at a 
high enough redshift to be able to achieve (7„ = 3 cm s^^ or better 
on each object individually, then we would require O ~ 2.1. Note 
however, that the ordinate of Fig.ll2ldoes not give the overall value 
of (Ji, for a i experiment using the best Nqso targets and setup O. 
The reason is of course that the total value of CTi, of such an exper- 
iment depends on how the total integration time is split up among 
the Nqso targets (see below). 

7.4 Target selection and observing strategies 

The question of how to select targets for a z experiment depends 
on what exactly one wishes to achieve. For example, we could sim- 
ply rank the candidate targets by their achievable CTi, for a fixed 
observational setup O, as suggested by Fig.[TT] and select the ob- 
jects with the smallest values. This would be the correct selection 
strategy if the goal were to obtain the smallest possible overall (7„ 
of the combined sample. Another possible goal might be to max- 
imise the significance with which z ^ can be detected. In this 
case we would select the targets with the largest values of |«|/cri,, 
which would give much larger weight to the high redshift objects 
and would practically deselect all objects at zqso < 2.5 where the 
redshift drift is too small to be detected. Yet another option would 
be to maximise the sensitivity to Qa (at a fixed, given value of JIm) 
by selecting the objects with the largest values of It is 

clear that the choice of selection method will also depend on the 
possible adoption of priors from other cosmological observations. 

Once the exact goal and hence the target selection strategy has 
been defined we must also choose the number of targets to include 



in the experiment. Including more than just the 'best' target (ac- 
cording to a given selection strategy) is desirable for several prac- 
tical reasons, including the ability to observe for as large a fraction 
of the total available telescope time as possible (to ease schedul- 
ing), and to be able to identify any potential directional systematics 
caused by the atmosphere, the telescope, the instrument, the cal- 
ibration procedures or the transformation of redshifts to the cos- 
mological reference frame. Obviously though, the more objects are 
included into the experiment the worse the final result will be (in 
the absence of systematics) because observing time will have to be 
redistributed from the 'best' object to the less suited ones. 

Finally, we must also decide how to divide the total avail- 
able integration time among the objects included in the experiment. 
Again, there are several possibilities, including: the time could be 
split equally among the targets or in such a way as to ensure equal 
S/N, equal (7„ or an equal value of the selection parameter. (Note 
that in practice it is likely that operational and scheduling con- 
straints would limit the available choices.) 

From the above it is clear that there are a large number of dif- 
ferent possibilities of implementing a z experiment. In particular, 
the decision on what precisely the objective of such an experiment 
should be does not appear to be straightforward. We will now focus 
on three alternatives which may be perceived to be representative 
of three different approaches. 

1. The simplest possible goal is to aim for the most precise 
2 measurement possible, i.e. the smallest overall (t„. This may be 
considered a 'pure experimentalist' approach where virtually no 
prior observational information or theoretical expectation is used 
to (mis-)guide the design of the experiment. 

2. Another approach is to emphasise the most basic and yet 
most unique (and perhaps most captivating) feature of the redshift 
drift experiment, which is being able to watch, literally and in real 
time, the Universe changing its rate of expansion. In this case the 
aim is to prove the existence of a dynamical redshift drift effect 
and hence to measure a non-zero value of z with the highest possi- 
ble significancelfl However, selecting targets according to this ap- 
proach requires the adoption of, and hence dependence on, a spe- 
cific model for z{z), and is therefore not free of prior assumptions, 
in contrast to our first example. 

3. As explained in the introduction, the discovery and the un- 
known physical source of the acceleration of the universal expan- 
sion provide a strong motivation for any observation seeking to de- 
termine H{z), and this is also true for the direct and dynamical 
method of measuring z. From this point of view the goal should 
be to place the strongest possible constraints on the acceleration. 
Translating this goal to a target selection strategy again requires a 
model of z{z), including the acceleration. 



7.5 Simulations of example measurements and constraints 
on model parameters 

We now proceed to implement each of the above approaches on 
the QSO sample shown in Fig.[TT]in order to illustrate what can 
be achieved with each of them. However, we stress that the above 



The value of this approach is perhaps best appreciated by picturing the 
PI of a future i experiment angrily muttering "And yet it moves!" under 
her/his breath as s/he leaves a seminar room after having tried to convince 
a sceptical audience that (i) an effect has indeed been observed and (ii) that 
it is not an instrumental artifact. 
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Figure 13. The colour image and the contoure show the final, overall value 
of (Ti, achieved by targeting the A^qso objects with the best individual Ci, 
values (cf. Fig. II U and by employing a given combination of telescope 
size, efficiency and total integration time. The total integration time is split 
equally among the targets. The contour levels are at cj;"' = 2, 3, 4 and 
5 cm s~^. 



approaches and the specific implementations below are only exam- 
ple strategies, and that many variants, refinements and completely 
different procedures are possible. 

We begin with the simplest case of selecting targets by the 
value of their achievable (t„, which can be read off Fig.[TT]for each 
object. Having chosen the best A^qso targets by this criterion we 
can compute the total, overall value of (t„ that can be achieved for 
this sample and for a given observational setup, O: 



E 



{zQSO,i,mi, /obj.i O) 



(28) 



where /obj.i is the fraction of the total integration time allocated 
to the ith object. In Fig. [T3] we plot as a function of Nqso 
and O, where we have chosen to distribute the time equally among 



all objects (i.e. /( 



obj.i 



N. 



QSO 



). Distributing the time to give equal 



S/N instead (which may be operationally desirable) leads to slightly 
better values of al°^ when only a few QSOs are targeted. The rea- 
son is of course that the object with the best a„ is only the seventh 
brightest target (cf. Fig. II U and so it is allocated relatively more 
time. However, for Nqso ^ 15 this advantage is lost and equal 
time allocation produces al°^ values that are smaller by -^10 per 
cent compared to equal S/N allocation. 

From Fig.[T3]we can see that an overall value of al°^ ~ 2- 
3 cm s~^ is well within reach of an ELT, even when tens of objects 
are targeted for the experiment. The steepness of the contours also 
indicate that cr*°' does not depend very sensitively on Nqso, at 
least for iYqso > 10. 

For further illustration we now perform a MC simulation of 
a redshift drift experiment using this target selection strategy. We 
(arbitrarily) choose Nqso ~ 20 and O = 2 so that the overall 
accuracy of the experiment is cr,I;°' = 2.34 cm s~^ (cf. Fig. |13l>. 
Furthermore we assume that the observations span a time interval 
of Ato = 20 yr and that the true z{z) is given by our standard 
cosmological model. The result is presented in Fig. [14] where we 
show as blue dots one realisation of the 'observed' velocity drifts 




Figure 14. The three sets of 'data' points show MC simulations of the red- 
shift drift experiment using the three different example implementations 
discussed in the text. In each case we have assumed an observational setup 
of O = 2 and we plot as 'data' points the 'observed' values and er- 
ror bars of the velocity drift ii, expected for a total experiment duration 
of Ato = 20 yr and for standard cosmological parameters (hro = 1, 
fijvi = 0.3, = 0.7). For a given QSO we use the centre of the Ly« 
forest as the redshift of the v measurement. Blue dots: selection by ay, 
Nqso = 20 (binned into four redshift bins), equal time allocation. Yellow 
squares: selection by \v\/cry, Nqso = 10 (in two redshift bins), equal 
time allocation. Brown triangles: selection by best combined constraint on 
SIai ^QSO = 2, optimal time distribution. The solid lines show the ex- 
pected redshift drift for different parameters as indicated, and h^o = 1- 
The grey shaded areas result from varying Ho by ±8 km Mpc~t 



and their errors along with the input model (red solid line). Since 
the selected 20 targets cover the redshift range 2.04 ^ zqso 5; 
3.91 quite homogeneously we have binned the measurements into 
four equally sized redshift bins. 

By construction these points represent the most precise mea- 
surement of i that is possible with a set of 20 QSOs (using equal 
time allocation) and O — 2. However, since many of the selected 
QSOs lie near the redshift where z — (for the assumed model) 
the redshift drift is only detected with an overall significance of 
5 — \v\/al°^ = 1.4, where \v\ is the weighted mean of the abso- 
lute values of the expected velocity drifts. 

This can be improved upon by choosing our second approach 
and selecting targets by the largest value of \v\/ai,. This quantity 
is a strongly increasing function of redshift and this selection strat- 
egy results in quite a different set of objects: the best 20 targets 
according to this criterion include only 3 of the objects previously 
selected. For O = 2, Ato = 20 yr and equal time allocation as 
before, the combined overall significance of this sample is 5 = 2.7 
(assuming our standard model of z is correct). The yellow squares 
in Fig.ll4lshow the result of a MC simulation using this setup and 
implementation, except that we use A^qso ~ 10 (in two redshift 
bins) which gives a slightly better significance of 5 = 3.1. To reach 
<S ^ 4 we need to further reduce Nqso to 3, or instead increase O 
to 3.4 or Ato to 25.4 yr. 

Finally, we turn to our third approach and the question of how 
to best select targets to constrain the acceleration of the expansion 
and what can be achieved in this respect with our sample of known 
QSOs. As in the previous case the answer will depend on what to 
expect for the expansion history and in particular for the acceler- 
ation. For the purpose of the following discussion we will again 
assume our standard cosmological model. 
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The simplest thing we can do to constrain the acceleration is to 
unambiguously detect its existence, i.e. to measure i > with the 
highest possible significance. This implies (i) that we need a v mea- 
surement at z < zo, where zq is defined by v(zo) — 0, and (ii) that 
target selection should proceed by the largest value of v/av, which 
indeed favours the lowest available redshifts. However, even if we 
use only the single best object by this criterion and assume a gener- 
ous Ato = 25 yr then a 2cr„ detection of v > would still require 
an unfeasible O — 9.5. The reason for not being able to do better 
is of course our inability to access the Lya forest at 2 ~ 0.7 where 
V is the largest. 

Let us analyse the situation more systematically by switch- 
ing to the parameter space of our cosmological model (which is 
3-dimensional since we will not assume spatial flatness). In this 
parameter space our goal of detecting the acceleration translates to 
proving that the deceleration parameter qo = — is < 0. 
However, in our model the acceleration is due to a cosmological 
constant, and so for simplicity we will instead pursue the slightly 
easier goal of proving the existence of a cosmological constant, i.e. 
of placing a positive lower limit on f^A after having marginalised 
over SIm and Hq. 

Consider the constraint of a single measurement of v — 
Av/Ato at some z (i.e. provided by some QSO), which happens 
to be precisely equal to its expected value. Obviously, a single data 
point cannot constrain all three parameters at the same time but 
only some 2-dimensional surface. Solving equations Q and l[5j for 
Qa we find that this surface is defined by: 



^A 



(1 + ^) 



(1 + 2)2 



z 



V 



(29) 



Since v is assumed to be equal to its expectation value this con- 
straint surface is entirely defined by 2 and it must include the point 
(JIm, i^A, ftro) ~ (0.3, 0.7, 1). Furthermore, the measurement er- 
ror (Ti, endows the surface with some small thickness. Hence we 
can see that each QSO in our sample [representing a (2, ct^,) pair] 
defines its own constraint surface. The question is whether there ex- 
ists some combination of some sub-set of these surfaces such that 
the projection of the combined constraint onto the f^A-axis does not 
include the region JIa ^ 0. If so, what is the largest significance 
level at which Q,a = can be excluded? 

To address this problem let us fix ^70 = 1 for now, so that 
a given constraint surface turns into a band in the fiA-f^M plane 
(cf. the grey shaded bands in Fig. 1151 ). Consider first the case 
Nqso = 1- For a fixed Hq, equation i29\ simply gives a linear 
relation between JIa and SIm which has a non-zero slope for all 
2; > 0. It follows that the projection of the constraint band onto the 
fiA-axis will be infinite unless we can place some kind of limit on 
JIm. Since the slope is positive we require a lower limit on Qm in 
order to obtain a lower limit on SIa. Rather uncontroversially we 
can impose JIm > 0. Applying this limit to equation ( 1291 ) we find 
that a single i) measurement can place a positive lower limit on f^A 
if the intercept of the fiA-f^M relation, 6, is significantly larger than 
0. The error on b is given by the width of the constraint band along 
the fiA-direction: 



Av cHq (1 + 2)2 



2 {l + z)H{z)_a^ 
1 Ato' 



(30) 



and the object giving the best lower limit on JIa is the one with the 
largest value of b/ Ab. However, we find that even the best object 
by this criterion does not deliver a significant lower limit on f^A. 
Hence we find that not only are we unable to detect v > with any 




Figure 15. Expected constraints in the Qa'^m plane from two measure- 
ments of the redshift drift at two different redshifts as indicated, assuming 
= 2 and Ato = 20 yr. The two targets and the spUt of the observing 
time among them were chosen to give the best lower Umit on Qa- The red 
and blue solid lines and the grey shaded bands show the individual con- 
straints provided by each of the two objects assuming a fixed hro = 1. The 
coloured ellipses show the joint 68 and 90 per cent confidence regions that 
result from combining the two measurements, marginalising over Ho using 
an external prior of Hq = (70 ± 8) km s~^ Mpc^^. The hashed region 
indicates the 95 per cent lower limit on Qa- The dashed contour shows 
the 68 per cent confidence region that results from using a flat prior on 
Ho and just marginalising over the range ^ Hq ^ 140 km s~^ Mpc^^ . 
Flat cosmologies and the boundary between current de- and acceleration are 
marked by solid black lines. The dark shaded region in the upper left comer 
designates the regime of 'bouncing universe' cosmologies which have no 
big bang in the past. 



significance but it is also impossible to prove D,a > with a single 
V measurement. 

However, this conclusion rests entirely on the weak lower limit 
on Hm that we used (flu > 0). The larger this limit the stronger 
also the lower limit on Qa- In other words, a given v measurement 
translates to a different constraint on Qa depending on what we 
assume for Qm ■ Hence, even if we cannot detect i; > with high 
significance it is nevertheless possible to place a significant lower 
limit on Qa- All we need is a stronger constraint on Qm- 

The point is that this constraint can be supplied by the 2 exper- 
iment itself, without having to resort to an external prior, by means 
of a second measurement at high redshift (as opposed to at 2 < 20). 
The idea is to combine the measurements from two objects (or two 
groups of objects) that are selected in different ways, such that their 
constraint bands in the SIa-S^m plane make as large an angle as pos- 
sible while also being as narrow as possible (cf. Fig.llSt. The first 
group is selected to have shallow constraint lines, implying low 
redshifts (2 < 20) and supplying a lower (but by itself insufficient) 
limit on Qa as discussed above, while the second group is selected 
to have steep constraint lines, implying high redshifts (cf. equa- 
tion |29) and supplying a lower limit on Qm- This is illustrated in 
Fig.fUJwhere the red and blue lines and the grey shaded areas show 
the constraint bands from two v measurements at low and high red- 
shifts as indicated. If the angle between them is large enough and 
if the individual a^s are small enough then the combination of the 
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two constraints (in a sense) results in an ellipse in the f^A-f^M 
plane whose lower edge excludes Q.a ~ 0. 

Hence we find that we need at least Nqso — 2. However, re- 
call that so far we have kept Ho fixed. Let us now reconsider the 
above as a function of Ho- First of all, we note that the constraint 
line resulting from the v measurement at z ~ zq is very insensitive 
to the value of Hq (because w « 0). Secondly, for the high redshift 
measurement the constraint line shifts upwards in Fig.[T5]as Ho in- 
creases. Hence the centre of the joint ellipse (where the constraint 
lines cross) moves from (Om,J^a) = (0.3,0.7) approximately 
along the low-z constraint line towards the lower left of Fig. [TS] 
However, at the same time the extent of the ellipse decreases be- 
cause a larger Ho predicts a larger redshift drift and hence results 
in stronger constraints on and for fixed measurement er- 
rors. The net result is that the lower edge of the joint ellipse and 
hence the obtainable lower limit on SIa are relatively insensitive to 
Ho- As we will see more quantitatively below, a third v measure- 
ment to constrain Ho is therefore not necessary. 

Thus we find that the strategy of selecting two groups of tar- 
gets as described above should provide the best possible lower limit 
on Qa, even after marginalising over Ho - However, when trying to 
express this approach in equations in order to proceed with the tar- 
get selection one realises that the selection of the two groups cannot 
in fact be separated, making an analytical procedure impractical. 
Hence we resort to the simplest case of only one object per group 
and simply try all possible combinations of two objects in our sam- 
ple. For any given combination we also determine empirically the 
optimal split of the total integration time among the two objects. In 
this way we find that the best lower limit on Qa is achieved for the 
two very bright objects in Fig.[TT]at zqso = 2.15 and 3.91, with 
the latter object receiving 0.3 of the total integration time. Mock 

V measurements for these objects are shown as brown triangles in 

In Fig.[T5]we present the expected constraints from these two 
objects in the JIa-^^m plane, assuming O = 2 and Ato = 20 yr. 
The red and blue solid lines show the individual constraint lines for 
a fixed hro = 1 (equation |29t and the grey shaded areas around 
them are the error bands that correspond to the respective values of 
a„ for each object (equation [30l ). Note that the optimal time split 
(in the sense of giving the best lower limit on Qa) is such that the 
two bands are of roughly equal width. 

The coloured ellipses in Fig. [TS] show the joint 68 and 90 per 
cent confidence regions that result from the combination of the two 

V measurements. Here we have assumed that an external constraint 
place s la limits of ±8 km s~^ Mpc~^ on Ho ( Freedman et al.l 
l200lh before marginalising over this parameter. By comparing the 
error bars on the brown triangles in Fig. [14] with the grey shaded 
region around the red line we can see that this uncertainty in Ho is 
small relative to the errors on v. Hence the ellipses in Fig. 1 151 are 
almost those that would be obtained for a fixed Ho- 

Projecting the joint constraints onto the vertical axis we find 
that the 95 per cent lower limit on the cosmological constant is 
JIa > 0.16 (indicated by the hashed region in Fig. I15t. while 
SIa = is excluded at the 98.2 per cent confidence level. Note 
that we can also exclude go = at a similar level of 97.6 per 
cent. These numbers do not depend sensitively on the adopted prior 
on Ho, for reasons explained above. Doubling the error on Ho to 
±16 km s^^ Mpc^^ has almost no effect at all. Even if we use 
a flat prior and simply marginalise over the range ^ Ho ^ 
140 km s~^ Mpc~^ we can still reject f^A = at 98.1 per cent 
confidence, although the 95 per cent lower limit now reduces to 
0.08. For a shorter experiment duration of Ato ~ 15 yr (and the 
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Figure 16. Expected constraints in the Qa-Ho plane for flat cosmologies 
assuming O = 2 and Ato = 20 yr and using the same two redshift drift 
measurements as in Fig. [15] The red and blue solid lines and the grey shaded 
bands show the individual constraints provided by each of the two objects. 
The coloured ellipses show the joint 68 and 90 per cent confidence regions 
that result from combining the two measurements. 



flat Ho prior) we can still exclude f^A = at 93.5 per cent confi- 
dence. 

Thus we conclude that, assuming general relativity to be valid, 
an ELT will enable us to confirm the existence of a cosmologi- 
cal constant using purely dynamical evidence - by measuring the 
redshift drift in QSO absorption lines over a period of ~20 yr - 
without assuming flatness, using any other external cosmological 
constraints or making any other astrophysical assumptions what- 
soever. 

Combining the constraints of a i experiment with those from 
other cosmological observations would of course lead to signifi- 
cantly stronger results. A detailed examination of this issue is, how- 
ever, beyond the scope of this paper. We will let it suffice to make 
two general points: (i) the z constraints will be unique and there- 
fore new in the sense that they cannot be provided by any other 
cosmological observation, i probes H{z) in a redshift range that is 
inaccessible to other methods (such as SNIa, weak lensing, BAO) 
and this will remain true even in the ELT era. In terms of charac- 
terising H{z), z measurements at z > 2 will nicely complement 
the data at z < 2 provided by SNIa surveys. In the fiA-f^M plane 
these datasets offer similar but nevertheless distinct constraints, 
(ii) The datasets that are most complementary to the redshift drift 
in the SIa-^m plane are those that constrain the geometry of the 
Universe, such as the fluctuation power spectrum of the CMB, 
which can break the remaining degeneracy in the i constraints. In 
fact, the WMAP 3 -y ear data constrain the quantity SI a + 0.72 JIm 
JSpergel et alj2007n , so that the CMB degeneracy line is almost ex- 
actly orthogonal to the i constraints shown in Fig.[T5] Combining 
these constraints would lead to individual 2a errors in Qa and Qm 
of ~0.08 and ~0.06, respectively. 

Finally, let us consider models in which the Universe is ex- 
actly flat (fift = 0). The parameter space of these models is 2- 
dimensional and a single v measurement results in a constraint line 
given by 
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In Fig.[T6]we show the constraints in the i^A-Ho plane offered by 
the two V measurements considered above (O = 2, Ato = 20 yr). 
Clearly, Hq is not constrained in any useful way and both substan- 
tially higher and lower redshift measurements of similar or better 
quality would be required to do so. On the other hand, JIa is con- 
strained quite well. Marginalising over Ho we find a 2a range of 
[0.42, 0.74], and Q,a = is excluded at 98.7 per cent confidence. 
For a shorter experiment duration of Ato = 15 yr we can still ex- 
clude Qa = at 94.7 per cent confidence. Note that these latter 
values are not very much higher than in the general case above. 
The difference is of course that in the flat case the best fit along 
the JIa = axis is obtained for an unfeasibly low Hq value of 
hro ~ 0.1, while in the general case this value is 1. Hence, for 
flat cosmologies any external Hq prior has a strong impact on the 
lower limit one can place on Qa, in contrast to the general case. 
For example, applying the relatively weak prior that Hq is known 
to within ±16 km s~^ Mpc^^ results in the rejection of Q,a = at 
> 99.99 per cent confidence for any Ato > 9 yr or, alternatively, 
for Ato = 15 yr and any O > 0.5. 



7.6 Constraints on alternative cosmological models 

So far we have only considered the standard cosmological model, 
i.e. general relativity with a stress-energy tensor dominated by cold 
(dark) matter and a cosmological constant (ACDM). Needless to 
say, industrious theorists have proposed numerous alternative ex- 
planations for the observed acceleration of the universal expansion. 
Broadly speaking, these models either suggest alternative forms of 
dark energy or modify the theory of gravity. In Appendix|B]we give 
a brief overview of two examples for each of these two classes of 
models. These are quintessence, generalised Chaplygin gas, Car- 
dassian expansion and the DGP brane-world model. Despite fun- 
damental differences, some generic features are shared by all four. 
Obviously, they all 'predict' late-time acceleration. They also agree 
that the early evolution of the Universe is indistinguishable from 
matter-dominated expansion in general relativity (so as not to spoil 
the successes of the standard model in explaining early structure 
formation). Furthermore, all of them have at least as many free pa- 
rameters as ACDM. In fact, three of the four models considered 
here include ACDM phenomenologically as a special case, i.e. as 
a sub-space of the models' parameter space. With these properties 
in mind, and given that ACDM can successfully reproduce SNIa, 
CMB and BAO data, it is perhaps not too surprising that the alter- 
native models can also be adjusted to fit these datasets, resulting in 
best-fi t parameter values consistent with ACDM (e.g. iDavis et aU 
l2007l) . 

Given this state of affairs we do not think it justified to study 
these alternative models, in this paper, to the same level of detail as 
we did for ACDM in the previous section. The name of the game 
would be to ask, for each model, how much of the non-ACDM-like 
parameter space could be excluded by a redshift drift experiment. 
This would best be done by considering the i constraints jointly 
with those from other relevant datasets. Instead, we w ill do some- 
thing m uch simpler here. As already mentioned above, lo^avis et al.l 
(120071) have fit all four of our alternative models to current SNIa, 




CMB and BAO data. The idea is now to ask, for each model, which 
values of i are still 'allowed' giv en the const raints already placed 
on the model's parameters by the lDavis et alj analysis. Comparing 



Figure 17. Comparison of the predictions for the redshift drift by non- 
standard cosmological modefs. For simpHcity, we impose flatness and fix 
hro = 1 in all cases, (a) At a given redshift the vertically hashed region 
shows the range of z values predicted by quintessence, allowing the model's 
three free parameters to vary within the join t 95.4 per cent co nfidence re- 
gion imposed by SNIa. CMB and BAO data ("P avis et al]|2007 l). The green 
and red solid lines show the predictions for those parameter combinations 
that produce the most extreme allowed z values at 2 = 0.3 and 4, respec- 
tively. (uiQ, Wa) = (— li 0) corresponds to flat ACDM. (b-d) Same as (a) 
for the other models as indicated. For the Chaplygin gas and Cardassian 
models a = and (g, n) = (1,0) correspond to flat ACDM, respec- 
tively. In each panel we also show, for comparison, the predictions of the 
flat ACDM model with 0.65 ^ f2A ^ 0.74 as the grey shaded region. We 
also reproduce the three sets of simulated measurements from Fig. 1 141 



these allowed ranges with the simulated z measurements of the pre- 
vious section will give a good impression of what can be achieved. 

In Fig.[T7]we plot this allowed z range as the vertically hashed 
region for all four models considered here along with the predic- 
tions for some specific parameter combinations as indicated. Here 
we have assumed Ato = 20 yr. In order to reduce the number 
of free parameters we are only considering flat models with Ho 
fixed at h7o — 1. For each model the plotted i range corresponds 
to the joint 95.4 per cent confidence region in the space spanned 
by the model's remaining free parameters. The grey shaded region 
in each panel shows the flat ACDM model (which has only one 
free parameter since we have fixed Qk and Ho), varying Qa in the 
range [0.65, 0.74]. This corresponds to the 95.4 confidence range 
obtained from the two redshift drift measurements discussed in the 
previous section. Note that the hashed and grey shaded regions are 
not 'centred' on each other. The reason is that we chose to centre 
the grey shaded region on our fiducial value of SIa = 0.7, whereas 
the ranges of the equivalent parameters in panels (a)-(c) are cen- 
tred on 0.73. For reference, we also reproduce, in each panel, the 
simulated data points from Fig. 1 141 

In panels (a)-(c) the hashed regions are somewhat larger than 
the grey shaded regions and the error bars on the 'data' points, in- 
dicating that a redshift drift experiment will be able to significantly 
reduce the sizes of the allowed regions in these models' parame- 
ter spaces. The specific parameter combinations that we chose for 
the solid lines are those that produce the most extreme allowed 
z values at z = 0.3 (green lines) and at z — 4 (red lines). For 
the quintessence and Cardassian models we can see that these two 
sets of lines have very different parameter values. Indeed, the two 
parameter combinations that produce the extreme z values at low 
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redshift predict very similar values at high redshift, and vice versa 
(although to a lesser extent). This again highlights the complemen- 
tarity between H{z) constraints at low and high redshifts. Notice, 
however, that the situation is different in panel (b). In the case of 
generalised Chaplygin gas one set of parameters produces the max- 
imum (or minimum) allowed i values at more or less all redshifts. 
The reason is that the regio n in this mo del's 2-dimensional pa- 
rameter space allowed by the lOavis et al.l analysis is in fact nearly 
1-dimensional. This is of course exactly true for the DGP model 
(panel d) which has only one free parameter to begin with. Notice, 
however, that in this case v{z) has a somewhat different slope. In 
contrast to the other models, the DGP prediction for i cannot be 
made arbitrarily close to the ACDM prediction and from Fig. [TT] 
we can see that the differences are of roughly the same size as the 
expected i measurement errors. 

We defer a detailed investigation of these issues to a future 
analysis. Here we simply conclude that the increased redshift range 
of H{z) measurements afforded by a redshift drift experiment will 
significantly add to our capability of discriminating between alter- 
native cosmological models. 



8 CONCLUSIONS 

In this paper we have appraised the prospects of using an ELT 
to detect and characterise the cosmological redshift drift, i, by 
means of a ~two decade-long spectroscopic monitoring campaign 
of high redshift QSO absorption lines. This experiment would be 
unique among cosmological observations in that it directly probes 
the global dynamics of the Robertson- Walker metric, thereby offer- 
ing a novel and truly independent path to the expansion history of 
the Universe. We summarise our results a s follows: 

1 . We confirm the earlier proposition ( lLoebll99^ that the Lya 
forest is the most promising target for a redshift drift experiment. 
We have explicitly examined the peculiar motions of the gas re- 
sponsible for the Lya forest in a hydrodynamical simulation and 
found their effects negligible. We have also assessed the effects of 
variations of the physical properties of the absorbing gas over the 
timescale of the experiment and found them similarly insignificant. 

2. Under the assumption of photon-noise limited observations 
the accuracy, Ov , to which a putative radial velocity shift between 
two spectra of the same object can be determined depends chiefly 
on the number and sharpness of the relevant spectral features and 
the spectra's S/N. The number density of the absorption lines that 
constitute the Lya forest depends strongly on redshift. Using exten- 
sive Monte Carlo simulations we have hence derived a quantitative 
relation between the (t„ of a pair of Lya forest spectra on the one 
hand, and the spectra's S/N and the background QSO's redshift on 
the other hand (equation [TS}. 

3. Apart from the Lya forest, all QSO spectra also display 
higher order H I lines and a variety of metal absorption lines. Using 
the H I Ly/3 lines in addition to the Lya lines, as well as all available 
metal lines within and to the red of the Lya forest significantly 
improves a given spectrum's sensitivity to radial velocity shifts by 
a factor of 0.67 (equation [T6}. 

4. In practice it will not be possible to measure the redshift 
drift by simply comparing just two observations of some set of ob- 
jects obtained at two well-defined epochs separated by some Ato- 
In reality, the total observing time required for the experiment will 
be spread over many nights (= epochs) whose distribution within 
Afo (now defined as the interval between the first and last obser- 
vations) will be subject to many practical constraints. In Section|6] 



we have calculated the relevant 'form factor' by which (t„ increases 
(or, equivalently, by which the effective Ato decreases) as a result 
of a given observing time distribution. For realistic distributions we 
find factors of ~ 1.1 to 1.7. 

6. The detailed implementation of a redshift drift experiment 
will depend on the exact objectives one hopes to achieve. We have 
elaborated three different possible top-level goals and we have 
worked out the specific target selection strategies required to tailor 
a i experiment to each of these. Applying these selection schemes 
to a comprehensive list of high redshift QSOs already known today 
we use our final scaling relation to predict an ELT's capability 
to measure i. For example, selecting those 20 targets that deliver 
the greatest measurement accuracy we find that a 42-m ELT could 
achieve a total, overall accuracy of 2.34 cm s~^ in 4000 h of ob- 
serving time (Fig.ll3t. If we additionally assume that our standard 
cosmological model correctly predicts the size of the redshift drift 
then we will be able to prove its existence (i.e. measure a non-zero 
value) at a significance of 3.1(t in Afo = 20 yr (cf. Fig.ll4t. 

Unfortunately however, due to the fact that the Lya forest is 
only accessible from the ground for z > 1.7 it will not be pos- 
sible to directly measure i > (the hallmark of past accelerated 
expansion) from any one of the QSOs in our sample with any sig- 
nificance, unless the transition from decelerated to accelerated ex- 
pansion occurred significantly earlier than predicted by our stan- 
dard model. However, even if the standard model's prediction for 
2 is correct an ELT will nevertheless be able to provide unequivo- 
cal proof of the existence of past acceleration by obtaining an upper 
limit on z at high redshift and a lower limit at 2 ~ 2 where the tran- 
sition from i < to > is expected to occur (cf. Fig.ll4t. For the 
above experiment parameters these limits will be strong enough to 
prove that 2 must become positive below 2 ~ 2. We stress that the 
validity of this proof only requires that gravity can be described by 
a metric theory and that the Universe is homogeneous and isotropic 
on large scales. It is entirely independent of any specific theory 
of gravity (such as general relativity) and does not require any as- 
sumptions regarding the spatial geometry of the Universe. It is also 
independent of other cosmological observations and does not rely 
on any other astrophysical assumptions whatsoever. Hence, a i ex- 
periment carried out with an ELT will arguably provide the most 
direct evidence of acceleration possible. 

In the context of general relativity dominated by cold dark 
matter and a cosmological constant [(nM,f2A) ~ (0.3,0.7)] a 
2 experiment over Ato — 20 yr using 4000 h of observing time 
on a 42-m ELT will be able to exclude Q,a = with 98.1 per 
cent confidence, providing a 2a lower limit of 0.08 (cf. Fig. 1151). 
Additionally assuming spatial flatness results in a 2a range of 
0.42 < < 0.74 (Fig.[T6ll. 

7. Finally, we point out that a i experiment constrains H{z) 
over a redshift range that is inaccessible to other methods. Redshift 
drift measurements at 2 > 2 would therefore complement current 
and future SNIa and BAO surveys. Adding 2 to the suite of mod- 
ern cosmological tests would thus provide additional leverage in 
constraining the parameters of both ACDM and alternative cosmo- 
logical models (cf. Fig.llTt. 

We stress that all the results summarised above were derived 
under the assumption that the uncertainty on radial velocity shift 
measurements will be dominated by photon noise. The next step 
is to understand to what extent this assumption can be realised 
in practice. To this end we need to (i) draw up a list of all po- 
tential sources of error, both astrophysical and instrumental in na- 
ture, (ii) quantitatively assess their impact on a 2 measurement, and 
(iii) where necessary and possible identify appropriate hardware 
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solutions, calibration and/or analysis techniques that will reduce 
their impact below the level of the photon noise. Although we have 
already addressed some of the most obvious astrophysical sources 
of noise in Section [3TT] a realistic and comprehensive assessment 
of the total error budget is a non-trivial task - in particular because 
of point (iii) above - which we defer to a future paper. 

One of the potential error sources that will require careful at- 
tention is the fact that QSO continua are to some extent variable. 
The flux-difference method we described in Section |43] to extract 
the redshift drift signal from the data is clearly very sensitive to 
any errors in determining the change of the continuum. Similarly, 
it is also quite susceptible to any residual uncertainties in the sky 
subtraction, flat-fielding and stray light correction. Clearly, a flux- 
difference type of approach would not be the method of choice for 
a real-life i experiment - a conclusion we had already reached in 
Section l431 for an entirely different reason. However, other methods 
of extracting the redshift drift signal are likely to be less sensitive to 
these issues. As already suggested in Section|5]we believe that the 
most promising approach is to simultaneously model the absorp- 
tion in all spectra of a given object with Voigt profiles, including 
the z of each component (or a model thereof) as an additional free 
parameter. The drawback of this method is that a complex absorp- 
tion feature can usually be modelled in more than one way, none of 
which necessarily has to correspond to reality, and this could rep- 
resent an additional source of error for a i measurement. In any 
case, the question of how to deal with QSO continuum variability 
remains to be investigated. 

Another issue we have not yet studied is the effect of the evo- 
lution of structure along the line of sight to a distant QSO. Consider 
two photons whose redshift difference we wish to measure. Along 
their way to us both will encounter various gravitational potential 
wells. Any evolution of these potentials between the passage of the 
first and second photon will cause them to follow slightly different 
paths, perturbing their redshift difference. This effect will have to 
be studied with high time resolution simulations of cosmic struc- 
ture formation. 

The quest for photon-noise limited radial velocity shift mea- 
surements also places rather exacting requirements on the spectro- 
graph to be used for a redshift drift experiment, in particular with 
respect to its stability and wavelength calibration. However, an in- 
strument concept study investigatin g these and other i ssues has al- 
ready been carried out (CODEX: iPasquini et al.|[200^ . concluding 
that the technical challenges are demanding but not insurmount- 
able. One of the novel concepts identified by this study is the use of 
a 'laser frequency comb' for wavelength calibration. This system 
provides a series of uniformly spaced, narrow lines whose abso- 
lute positio ns are known a prior i with relative precision better than 
10"^^ (see lMurphv et al.|[2007l for details). Note that an absolute 
wavelength calibration of this precision would allow the combina- 
tion of data from different ELTs and would guarantee the legacy 
value of the data by enabling future astronomers to measure i over 
much longer timescales. 

In deriving our results we have also used the magnitudes and 
redshifts of the high-redshift QSO population as it is known to- 
day. Note that we have used all known QSOs, irrespective of their 
declination. In exchange, so to speak, we have refrained from spec- 
ulating on future discoveries of bright, high-redshift QSOs. It is 
clear that existing QSO catalogues do not yet cover the entire sky 
uniformly with the same completeness, even at the bright magni- 
tudes we are interested in here. Indeed, only one third of the ob- 
jects shown in Fig. [TT] lie in the south. Furthermore, since colour 
selection is one of the most popular methods of QSO candidate 



identification it is likely that our current sample is incomplete in 
the redshift range 2.4 < z < 3 where the colour s of QSOs are 
very similar to those of main sequence stars (e.g. [Richards et al.l 
[2002). Hence it is not unreasonable to believe that at least a few 
bright QSOs still remain to be discovered, especially in the south. 
Considering that a z experiment would benefit greatly from any 
new discoveries we propose that a complete, systematic survey for 
bright, high-redshift QSOs be carried out once the site for an ELT 
has been chosen. 

In addition, we propose that all potential target QSOs should 
be monitored photometrically. Not only would this remove the un- 
certainties of our present photometry (note that we had to rely on 
photographic photometry from the SuperCOSMOS Sky Survey for 
most of the brightest objects in Fig. II It but it would also allow us to 
account for QSO variability in the target selection. If the monitoring 
frequency were high enough (as can be expected from possible fu- 
ture facilities like the Large Synoptic Survey Telescope) then these 
data could also be used to time the spectroscopic observations of a 
given target to coincide with a peak in the object's lightcurve. 

We close by observing that any QSO spectra collected as part 
of a i experiment could also be used to address numerous other 
important issues, including the variability of the fine-structure con- 
stant and the primordial abundance of deuterium. Given the sub- 
stantial observing time requirements for a z experiment it is clear 
that these synergies will play an important role in any realistic at- 
tempt to establish a future ELT redshift drift campaign. 
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APPENDIX A: REDSHIFT DRIFT INCLUDING 
PECULIAR MOTION 

Equation Q for the redshift drift is only correct if both the observer 
and the source are comoving. In reality, both will have peculiar ve- 
locities and accelerations due to local gravitational potential wells. 
The purpose of this Appendix is to provide an exact calculation of 
the effects of peculiar motion on the observed redshift drift. 

Our own peculiar motion can be broken down into several 
components, such as the Earth's motion with respect to the solar 
system's barycentre, the solar system's rotation within the Galaxy, 
and the Galaxy's movement with respect to the Hubble flow. All of 
these components are or will soon be known with high prec ision so 
that t he observations can be corrected accordingly (Pasqui ni et al.l 
l2006l) . Hence we will not consider the observer's peculiar motion 
any further. 

Por a given object at large distance it is in general not possible 
to determine its peculiar motion or to correct its observed redshift 
for it. However, the peculiar motions of widely separated objects 
are expected to be uncorrelated. Hence, when using several such 
objects to measure i, their peculiar motions are only expected to in- 
troduce a random noise component, but will not cause any system- 
atic bias. To assess the relevance of peculiar motions for any given 
implementation of a i experiment one must compare the noise in- 
duced by the peculiar motions of the experiment's targets with other 
sources of error (ideally only photon noise). 

We now derive the expression for the observed redshift drift 
in the presence of peculiar velocity and acceleration. Consider a 
photon emitted at time ton by a source moving with velocity v 
with respect to its local Hubble flow (object A in Rig. lAlt . and 
received by a comoving observer at time tobs (object B in Rig. lAlt . 
Clearly, the observed redshift of this photon travelling from A to B 
is given by 

1 + -2obs(tobs, icm) = [1 + 2D(tcm)] [1 + 2(tobs,iom)]. (Al) 

Here, z is the cosmological redshift as before, and zu is the (special 
relativistic) Doppler shift of the emitter: 

l + zu{tcm) = ^ ^ " - =7(l + cosg-), (A2) 

where 9 is the angle between the emitter's peculiar velocity and 
the outward line of sight (LOS) from the observer (see Rig. lAlt . 
The redshift drift is observed as the redshift difference between the 
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Figure Al. Cartoon illustrating various photons, their emitters and ob- 
servers, as well as their times of emission and observation (see text for 
details). We also show various angles that define the peculiar motion of 
object A: the angle d in panel (a) is the angle between the peculiar velocity 
V and the LOS, while ri is the angle between v and the peculiar acceleration 
projected into the plane of the drawing, ap . The acceleration a itself makes 
an angle C, with this plane (not shown). In panel (b) v has acquired a com- 
ponent perpendicular to the plane of the drawing (making an angle ^), and 
so we only show its projected component. Dp, as well as its angle dp with 
the LOS. Compared to panel (a) v has shifted by an angle tp in the plane of 
the drawing while the direction of the LOS has changed by the angle <f>. 



photons labelled as 1 and 2 in Fig. lAll and has two terms: 



^obs (^ob 



! ''cm/ 



''obs 



*obs 



dtobs 



(1 + ^) 



Az 

dtobt 



(A3) 



We can therefore identify the following, separable effects of pecu- 
liar motions on the process of putting data points on Fig. [2] (i) The 
Doppler shift in equation JAU introduces an uncertainty on where 
to place a given object along the a;-axis of Fig. |2] Uncertainties in 
the placement along the y-axis are caused by (ii) peculiar accelera- 
tion, (iii) the Doppler shift in the second term of equation l lA3t and 
(iv) the change in the comoving coordinates of the source. This last 
point can in turn be separated into two distinct effects: the decrease 
of 6 (due to parallax) and the change of the distance between the 
observer and the source. The latter implies that the Az/Atohs ap- 
pearing in equation ( |A3t is not the same as that of equation 
which is only valid for a constant x- 

This is illustrated by photons 3 and 4 in Fig. lAll Photon 3 is 
emitted at tom by the comoving object C that is at the same loca- 
tion as the moving object A at time t\^, while photon 4 is emitted 
by comoving object D just as A is passing by at time ti„^. Photons 
3 and 4 are received by B at the same times as photons 1 and 2, 
respectively (i.e. t],^^ and tobs)- Clearly, the redshift difference be- 
tween these two photons is the cosmological shift we are seeking: 



Az = 2:D^B(tobs,*cm) " ^C^B (^obs , ^Li) 



(A4) 



To calculate Az we make use of photon 5 in Fig. lAll Like photon 3 
this photon is sent from C to B, although at a different time. Since 
C and B are comoving objects, the redshift difference of photons 3 
and 5 is given by equation We have: 



Az 



ZD- 



*B(iobs, 



2C->B(tobs, icm + Atcm) 



+ 2C^B(tobs, ^cm + Ate 



ZC-.B 



(^obs 1 ^cm ) 



dz 



Az\, 



Atob 



= (1 + 2) H(tom) Atcm 

+ [(1 + Z) i/(iobs) - H{Um)] AUbs 
= (1 + z) H{tobs) Atobs 

- (1 + 2) //(torn) (Ate - Afem) 

= (1 + 2) H{tobs) (tlbs — ilbs) 

- (1 + 2) H{tcin) (iom - icm), 



(A5) 



where Atom is the time required by photon 4 to cover the comoving 
radial coordinate difference, Ax, between C and D as seen from B. 
Ate is the time difference between the emission of photons 3 and 
5 by C, so that Atobs ~ (1 + 2) Ate - The last equality of equation 
dASI l holds by construction. Using 



cAtcm = a(tcm)Ax = COS 9v (ti^ — tl^) 



one obtains 



-5 1£2^ = [(1 + ^) (l + cos0^)]-i 



■'cm "cm 



(A6) 



(A7) 



and hence 

= {1 + z)H{tobs)- {I + cos 6 -y ^ H{t, 



{l+z)H{tobs) 



H{tcm)- 



(A8) 



1 +2D 

Turning now to the first term of equation iA3l the change in 
the Doppler shift is given by: 

d2D = ^ {tL - tL) (A9) 
so that 

— — = __ [(1 + 2) (1 + cose-)] \ (AlO) 

where Azo / dtom is the rate of change of the Doppler shift as mea- 
sured by a local comoving observer along the LOS to B. Using a 
prime to denote differentiation with respect to tom, we have 



2d = — (1 + 2d) , 

c \ 1 + cos 



cos 9 iv 
^ + 7 - 



7 sin 6* - B' . (All) 

c 



In general, the peculiar acceleration, a, can have any orienta- 
tion with respect to v. We denote the direction of a by the angles 
^ and r\. The former is the angle between a and the plane spanned 
by V and the observer (i.e. the plane of Fig. lAlt . The latter is the 
angle between Op, the projection of a into this plane, and v [see 
Fig. lAlf a)!. Hence we have 



V = cos 77 Op — cos 77 cosi^ a. 



(A 12) 



The angle Q changes with time because both the viewing angle 
and the direction of v change. In particular, v may acquire a com- 
ponent that is perpendicular to the plane of Fig. lAll so that it makes 
some angle ^ with this plane. Writing for the angle between the 
LOS and the projected velocity (see Fig. lAll l, we have 



COsO{tcYn) = cos 6lp (tern) COS^(tcm) 

and hence 



(A13) 



sin^ 6' = cosf sin6'p 6'p + cos6'p sin^ (A14) 

However, since this equation should be evaluated at tl^ and since 
^(^em) = and 0{tl^^) = 6p{tl^) by construction, we can simply 
ignore the perpendicular component and write 0'{tl^) = 6'p{tl^) 
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[unless 6{tl^) — 0, in which case the second term of equation 
dAllb is zero anyway]. 

The change in Op is given by the angles and between the 
dashed and solid lines in Fig. lAU b). Clearly, <f) is the parallax angle 
for a baseline of s'mOv {tl^ — tem)> the transverse distance 
between C and D, while denotes the change of u's direction in 
the plane of the drawing. Hence we obtain 



sin rj cos — 

V 



• sin r] cos — 

V 

1+2 



rfpx(tem,x) 



(A15) 



rfpx(tobs, X) 

dpx is the parallax distance betwee n two objects separated by a c o- 
moving coordinate distance x (e.g. IPeeblesll 1 993l : IPeacocldl 1 999h : 

dpx(t,x)=a(i)f^^^ 



where we have used the notation 

{sin X k = 

X k = 

sinh X k = 



+ 1 


-1 



(A 16) 



(A17) 



and its cosine equivalent H(x) ~ -^/l — fcS2(x). 

Collating all of the above we finally find for the observed red- 
shift drift: 



dtobs 



7 COS ( COS rj 



cos f 



1 + cos ( 



+ 7 



2V 



a sm a 

7 cos C sm 77 — 

' ^ ' c l + cos( 



l + z sin- e ^ 



+ 7 

dpx 1 + cos f 

+ (1 + Zohs) H{tohs) 
(1 + Zobs) Ho 

— 7 < H[z) — cos 1^ - 



cos{9 + rf) 



1 + cos t 



+ cos 77 7 



2^^ 



+ 



1 + Zobs 

dnxfz) C 



1 + Zobs) -ffo 
1 + 2obs V[ 



H{z) 



+ 



dp^{z) 



where we have used the notation 

sin 6 V 



Vt 



i + coseiii 

c 

for the emitter's apparent transverse velocity and 

flr — COS cos{6 + rj) a 



(A18) 



(A19) 



(A20) 



for the magnitude of a's component along the LOS. Comparing 
with equation we can see that the main difference is the modifi- 
cation of H{z) with a term proportional to the change in the radial 
velocity. The second order term stems from taking the derivative of 
7 in equation dAllt . while the last term is due to the parallax effect. 
Clearly, equation JAlSt reduces to equation Q for v — a = 0. 
Note also that ( is entirely degenerate with a and so we will set 
^ = from now on. 

The formulae for many special cases such as a = (con- 
stant motion), 9 = (radial motion) and 6 = n/2 (transverse 
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Figure A2. Top panel: the solid line shows the cosmological redshift drift, 
z{z), expected for comoving sources (cf. Fig.|2j. The three grids show how 
an object at 2 = 2, 3 and 4, respectively, moves in the observed Zobs'^obs 
plane as a function of v and a, assuming 8 = rj = = 0. For each grid, 
the vertical lines are lines of constant v, while the tilted horizontal lines are 
lines of constant a. The values of v and a are indicated along the middle grid 
in units of 10^ km s~^ and 10~* cm s~^. The other panels show the same 
for different values of 8 and rj as indicated (and = 0). For comparison, 
the grey shaded areas in each panel are the same as in the top panel. The 
insets show a blow-up of the region around the 2 = 3 point. 
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motion) are immediately apparent from equation l lAlSb . Another 
special case is that of a freely moving source. In the absence of lo- 
cal gravitational potential wells we have v oc a^^ (and C, — -q = G) 
so that a = v' — —vH{tcm)- Neglecting terms of order {^)'^ and 
higher, we obtain 



^obs 
1 + 2:0b 



Ho 



H{z) 



(1 + ^) + 



dpx(z) 



(A21) 



which shows that in this case purely radial peculiar motion pre- 
serves the form of equat ion (|4). Setting = 0, we recover the 
formula o f lTeubeil l fTgS^ : 



Zohs 



1 + Zobs 



Ho[l - Vl + ^MZ (1 + ^)] + 



dpx(2:) C 



(A22) 



Given an object's cosmological redshift z, and the parameters 
of its peculiar motion (v, 9) and (a, rj) we can now predict its ob- 
served position in a iobs vs. Zobs plot, i.e. the observed equivalent 
of Fig.|2] This is shown in Fig. lA2l for three objects at 2 = 2, 3 and 
4, respectively. The solid line in each panel shows the cosmological 
z{z) that one expects in the absence of peculiar motions. In each 
panel we fix 9 and rj at some value and use small grids to show how 
the three objects are perturbed away from their 'correct' positions 
(marked as dots) as a function of v and a. The grids cover the range 
-1 to +1 X 10* km s"^ in v and -3 to +3 x 10"** cm s"^ in a. To 
relate these scales to one another we point out that a constant ac- 
celeration of 1.5 X 10~* cm s~^ would result in a peculiar velocity 
of 10'' km s-i by 2 = 3. 

Obviously, the largest possible effects are obtained when 9 — 
7] = (top panel), and we reproduce this case as grey shaded areas 
in all other panels for comparison. The extent of the grids in the v 
and a directions can be scaled independently and almost arbitrarily 
by appropriately adjusting 9 and ry. However, when 9 + rj ^ n /2 
purely relativistic and/or second order effects take over. These are a 
factor of ~ 10-100 smaller in magnitude and can significantly dis- 
tort the grids. Note that since the grids are generally not symmetric, 
it is in principle possible (contrary to our initial expectation) that 
peculiar motions could bias a 2(2) measurement, even when using 
a large sample of objects whose peculiar motions are uncorrelated. 

Given typical values for the peculiar velocities and accelera- 
tions of a particular set of objects (e.g. Figs.[3]and|4l(, we can essen- 
tially read off Fig. lA2l the scale of the noise on a 2(2) measurement 
introduced by the peculiar motions. This must be compared to the 
photon noise on 2(2) for an individual object. In the case of the 
Lya forest, as discussed in this paper, the photon noise error on 2 
for an individual absorption line is about four orders of magnitude 
larger than the error induced by peculiar motions of the absorbing 
gas. We conclude that peculiar motions will only be relevant for a 2 
experiment if it is based on very precise measurements on a small 
number of objects with large peculiar velocities and accelerations. 



leave the general relativistic field equation intact and (ii) models 
that modify the theory of gravity. 

Recall that the validity of the Robertson- Walker metric was 
the only assumption we needed to derive equation Q for the red- 
shift drift (cf. Section|2l(. None of the models considered in this Ap- 
pendix violate this assumption, and so all we need to do in order to 
determine a given model's pr ediction for 2 is to specify the model's 
Friedman equation (see also ISzvdlowski. Kurek & Krawiedl2006l 
for an overview). 



Bl Alternative dark energy models 

The general relativistic Friedman equation as given in equation ^ 
assumes a specific form for each mass-energy component's equa- 
tion of state, namely pi = WiC? pi, where pi and pi are the compo- 
nent's pressure and density, respectively. Without this assumption 
the Friedman equation reads: 



H{z) 



Ho 



^+fi,(l+2)^ 



(Bl) 



where pm is the present-day value of pi. For an arbitrary equation 
of state, pi{pi), the evolution of pi is obtained by solving the gen- 
eral relativistic energy conservation equation: 



dpi 
da 



Pi + 



Pi{pi) 



(B2) 



Bl.l Quintessence 

Here, we take the term quintessence to simply mean a smooth, 
time varying dark energy component with equation of state px ~ 
wx(? px, where wx may also vary with time. In principle, wx and 
its evolution should be given by the scalar field equation of the un- 
derlying physical theory of this component. In practice, however, it 
is common to use a simple parameterisation of wx in ord er to be 
less model dependent. Here we will use the expression of iLindeJ 
(2003): 



Wx{z) = Wo + Wail —] =-"^0+ Wa—^ • 

V ao / 1+2 

In this case the Friedman equation is given by 

2 



(B3) 



H(z) 



Ho 



= Vlx(l + z) 



+ Qm{1 + zf + + zf 



(B4) 



where we have included a matter component with wm — 0. Note 
that for (wq, Wa) = (—1,0) the above reduces to the standard 
ACDM case. 



APPENDIX B: SUMMARY OF SELECTED 
ALTERNATIVE COSMOLOGICAL MODELS 

There exist a number of cosmological models that may offer an 
alternative to the standard ACDM scenario. The purpose of this 
Appendix is to provide a summary of the predictions for the red- 
shift drift of an incomplete subset of these alternative models. All 
of them were introduced specifically to explain the observed accel- 
eration of the universal expansion, and they can be separated into 
two distinct classes: (i) models that modify the stress-energy ten- 
sor of the Universe to include a new dark energy component, but 



B1.2 Generalised Chaplygin gas 

In the ACDM and quintessence models it is the different evolution 
of the densities of the two mass-energy components that causes the 
transition from a decelerating (matter-dominated) to an accelerat- 
ing (A/quintessence-dominated) Universe. In contrast, the gener- 
alised Chaplygin gas (GCG) model seeks to explain this transition 
with only a single background fluid. This leads to a rather exotic 
equation of state for GCG: 

p, = -Ad'p-", (B5) 
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where A and a are constants with A and a ^ —1. Clearly, 
for a — —1 GCG is identical to A and therefore the model is 
identical to the de Sitter Universe (containing A but no matter), 
while a — 1 correspond s to the standard Chaplygin gas model 
jKamenshchik et alJboOlT) . The above equation of state results in 
the following Friedman equation: 



H{z) 



Ho 



nc[B + (l-B)(l- 



\3(l+a)il/{l+c) 



(B6) 



where B = A/p];^". Note that we are neglecting the baryons 
which had been subsumed in f^M previously. In this case GCG 
is the only mass-energy component and we have flc — 1 — Qk- 
Inserting the solution for pc into the equation of state we obtain: 



Pc = 



-Bc^ 



B + {l-B){l + z] 



S(l + a) 



Pc = —Wciz) C pc 



(B7) 



From the above equations we can see clearly that for < B < 
1 the Hubble expansion and the GCG equation of state smoothly 
change from being matter-like at early times to being A-like at late 
times (possibly in the future), where the parameters B and a jointly 
control when this transition occurs and how long it takes. For the 
extreme values of B = lUc(O) = and 1 GCG behaves like pure 
matter and A, respectively, regardless of the value of a. Note that 
for Q = the GCG model reduces to the standard ACDM case. 



B2 Alternative gravity models 

B2.1 Cardassian models 

In CardassiaiQ models the Universe is assumed to be flat and to 
contain only matter but no dark energy. The acceleration of the ex- 
pansion is instead achieved by altering the Friedman equation. In 
the m odified polytropic Cardassian models the Friedman equation 
reads dWang et alj2003t) : 



H{z) 



Ho 



^M{i + zy[i + {Q.i^ -i){i + z) 



3g(n-l)i 



,(B8) 



where the parameters g > and n < 2/3 are constant. 
The case q = 1 corr esponds to the original Cardassian model 
jFreese & Lewis! l2002h which is, however, as far as the Hubble 
expansion is concerned, exactly equivalent to a flat quintessence 
model with wx = n — 1. Note that in Cardassian models spatial 
flatness is not equivalent to SIm = 1- 

At early times the Cardassian Friedman equation is dominated 
by the first term inside the brackets, resulting in normal, matter- 
dominated expansion. However, at late times the non-standard term 
dominates and gives rise to the acceleration. Such modifications 
to the Friedman equation may be motivated by self-interacting 
dark matter or by considering our observable Universe as a 3- 
dimensional brane embedded in a higher dimensional space-time. 



distances larger than some cross-over scale re- On scales smaller 
than Tc an observer on the brane thus measures normal (3-1-1)- 
dimensional gravity, while on larger scales gravity weakens, caus- 
ing accelerated expansion beyond re- If rc is chosen large enough 
then it will only enter the Hubble horizon at late times so that the 
early evolution of the Universe remains unaltered compared to stan- 
dard, matter-dominated general relativity, while at late times the 
leaking of gravity to the extra dimension accelerates the expansion 
without requiring any form of dark energy. 

Defining an rc-induced density parameter. 



(B9) 



we can write the Friedman equation of the DGP model as 
jPeffavet et alj2002h : 



\H{z)] 
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Ho 





!^r, + \/^T, + ^m{1 + Zf 



+ nk{i + zY. 

Note that in this model the curvature parameter is given by: 

fife = 1 - (y^^ + V^'-e + !^m) , 

SO that flatness implies 
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(BID 
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B2.2 DGP model 

Anoth er bra ne-world scenario is the iDvali. Gabadadze & PorratH 
(DGP: l2000l) model. In this (4-l-l)-dimensional model gravity can 
spread (or 'leak') to the flat, infinite extra spatial dimension at 



With dark energy being 'bad' <Whitel2007h and Cardassians being 'ugly', 
all we need now is a 'good' theory to complete the famous trio. 



